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1.0  INTRODUCTION 

This  report  summarizes  the  tasks  completed  and  results  obtained  during  the 
period  1  February  1990  to  31  January  1993  under  ONR  contract  number  N00014-90- 
C-0071.  The  overall  goal  of  this  research  has  been  to  contribute  to  a  better 
understanding  of  the  processes  which  cause  variations  in  the  radar  backscatter  from 
the  ocean  surface  on  scales  of  1  meter  to  1  kilometer,  in  order  to  enhance  the  utility 
of  radar  remote  sensing  techniques  for  oceanographic  purposes.  The  specific 
objectives  were  to  (1)  develop  a  statistical  model  for  the  ocean  surface  which  includes 
the  effects  of  all  the  relevant  processes  which  influence  the  small-scale  surface 
roughness,  (2)  combine  this  model  with  an  appropriate  electromagnetic  scattering 
model,  and  (3)  test  the  predictions  of  the  combined  model  using  empirical  data  from 
laboratory  and  field  experiments. 

During  the  first  year  of  this  project,  work  centered  on  the  theory  of  wave- 
current  interactions  using  the  wave  action  spectral  transport  equation.  This  task 
included  both  the  development  and  exercise  of  models  in  order  to  make  pre- 
experiment  predictions  for  the  High  Resolution  field  experiment.  The  work  resulted 
in  the  publication  of  a  paper  in  the  Journal  of  Geophysical  Research ,  a  copy  of  which 
is  included  as  Appendix  A  of  this  report.  Additional  efforts  under  this  task  which 
have  not  been  published  in  the  open  literature  are  described  in  section  2  of  this  report. 

Subsequent  modeling  efforts  included  consideration  of  the  radar  backscatter 
problem,  including  a  comparison  of  two  implementations  of  the  two-scale  or 
composite  model  and  an  exploration  of  a  new  scattering  model  for  breaking  waves 
based  on  the  Kirchhoff  approximation.  The  results  of  this  effort  are  described  in 
section  3  of  this  report. 

During  the  second  and  third  years  of  this  project,  model  verification  studies 
were  carried  out  using  data  collected  during  the  High  Resolution  pilot  experiment 
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which  took  place  in  September,  1991.  A  Doppler  radar  was  deployed  on  the  LADAS 
platform  during  this  experiment,  and  data  from  this  instrument  was  analyzed  and 
compared  with  other  measurements  made  during  the  experiment.  Synthetic  Aperture 
Radar  data  collected  by  the  ERIM/NAWC  P3  SAR  system  during  the  High  Resolution 
experiment  was  also  analyzed  during  the  final  year  of  this  project.  These  efforts  are 
described  in  section  4  of  this  report.  Finally,  a  summary  of  the  work  accomplished 
and  results  obtained  during  this  project  is  given  in  section  5. 
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2.0  WAVE  SPECTRUM  MODELING 

The  ocean  surface  may  be  described  statistically  in  terms  of  the  wave  height 
spectrum  S(k,  <J>)  where  k  is  the  wavenumber  and  4>  is  the  propagation  direction. 
Although  this  description  does  not  necessarily  imply  anything  about  the  dynamics  of 
the  surface,  the  assumption  is  commonly  made  that  the  Fourier  components  which 
make  up  the  spectrum  may  be  regarded  as  freely  propagating,  linear  surface  waves. 
Spatial  variations  in  the  spectrum  are  then  considered  to  be  due  to  interactions  of 
these  waves  with  variable  surface  currents,  surfactants,  the  atmosphere,  and  other 
waves. 

The  interaction  of  short  waves  with  surface  currents  was  investigated  during 
the  first  year  of  this  project  using  the  wave  action  equation,  and  the  results  were 
published  in  the  paper  reproduced  in  Appendix  A.  Subsequently,  further  work  was 
done  on  the  net  source  function  which  appears  on  the  right-hand  side  of  this  equation. 
Some  of  this  work  was  presented  orally  at  the  Spring  1991  AGU  meeting,  but  since  it 
has  not  been  published  in  printed  form  it  is  described  in  more  detail  in  the  following 
section. 

2.1.  NET  SOURCE  FUNCTION 

The  primary  energy  source  for  ocean  surface  waves  is  the  atmosphere. 
Although  the  mechanisms  responsible  for  the  transfer  of  energy  from  the  atmosphere 
are  not  entirely  understood,  they  may  be  considered  to  fall  into  two  categories.  The 
initial  formation  of  waves  from  a  calm  surface  is  thought  to  be  caused  by  turbulent 
atmospheric  pressure  fluctuations  which  move  across  the  water  surface  at  the  same 
speed  as  the  waves  (Phillips,  1957).  The  energy  input  due  to  this  mechanism  is 
denoted  here  by  the  function  a  (k,  <J>)  .  Once  waves  of  a  certain  amplitude  are 
present  a  feedback  mechanism  begins,  leading  to  an  exponential  wave  growth.  The 
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rate  of  energy  input  due  to  this  mechanism  is  therefore  proportional  to  the  spectral 
energy  density,  with  a  constant  of  proportionality  denoted  by  0  (k,  4>)  . 
Measurements  of  this  quantity  exhibit  a  large  amount  of  scatter,  but  the  expression 


ps=0.25 


Pa 


U 


COS  -1 


CO 


given  by  Snyder  et  al  (1981)  is  widely  accepted  for  wavenumbers  near  the  spectral 
peak,  while  the  expression 


PP=0 . 04 


u  cos 


for  |$-$J<u/2 


suggested  by  Plant  (1982)  is  generally  considered  to  be  valid  at  higher  wavenumbers. 
Here  U  is  the  wind  speed  at  5  meters  above  the  surface,  u,  is  the  friction  velocity 
(see  section  2.4  below),  c  is  the  phase  velocity  and  w  is  frequency  of  the  waves,  4>  is 
the  wave  propagation  direction,  $w  the  wind  direction,  and  pa  and  pw  are  the 
densities  of  air  and  water,  respectively. 

In  order  for  the  wave  system  to  reach  a  state  of  equilibrium,  there  must  exist 
energy  dissipation  mechanisms  which  balance  the  energy  input  from  the  atmosphere. 
One  such  mechanism  is  due  to  the  viscosity  of  the  water.  The  rate  of  energy 
dissipation  due  to  viscosity  is  given  by  4 \k2  times  the  spectral  energy  density,  where 
v  is  the  kinematic  viscosity  of  the  water.  However,  this  mechanism  is  capable  of 
producing  an  equilibrium  state  only  if  4vk2>p ,  which  occurs  only  at  very  high 
wavenumbers  or  very  low  wind  speeds.  Mathematically,  a  state  of  equilibrium  can  be 
produced  by  introducing  a  rate  of  dissipation  which  is  a  nonlinear  functional  of  the 
energy  spectrum.  The  simplest  such  term  would  be  proportional  to  the  spectral 
density  raised  to  some  power  n>i .  For  n-2  the  net  source  function  can  be  written 
as 
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=  a  +  ($-4vk2)  S -y(ak*S2 
at 

where  y  is  a  dimensionless  constant.  The  steady-state  solution  of  this  equation  yields 
'he  equilibrium  spectrum 

S  (Jc  <{?)  -  P-4vJc2  +  y/(P-4vic2)  z  +  4ay<A>k* 

°  '  2y<t*k 4 

Although  the  form  and  magnitude  of  the  Phillips  growth  term  is  uncertain,  it  is 
assumed  to  be  the  dominant  term  for  wavenumbers  near  the  spectral  peak.  In  this 
region,  the  equilibrium  spectrum  is  then  given  by  approximately 

It  is  further  assumed  that  in  this  region  the  equilibrium  spectrum  has  the  form  given 
by  Pierson  and  Moskowitz  (1964).  This  would  imply  that  the  Phillips  or  linear 
growth  term  can  be  written  as 

a=y<i)/c45|M(/c,  <J>) 

where  SPM(k,$)  is  the  Pierson-Moskowitz  spectrum  converted  into  wavenumber 
space,  i.e. 

Sm(k,to)  =O.OD81gzu-5e'0-74("o/")V(<|>) 

where  g  is  the  gravitational  acceleration,  o >a=g/U,  U  is  the  wind  speed  at  19.5 
meters  above  the  surface,  F(4>)  is  an  angular  distribution  function  such  as 
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and  (a  is  given  by  the  dispersion  relation 


oiz=gk+  —  Jc3 
P 


where  t  is  the  surface  tension  and  p  is  the  density  of  water. 

At  high  wavenumbers,  where  the  viscous  dissipation  term  dominates,  the 
equilibrium  spectrum  approaches 


S0(k,4>) 


a 

4  vie2 


which  falls  off  as  jc-10  for  k2>pgr/x .  The  range  of  wavenumbers  where  this 
approximation  holos  is  sometimes  referred  to  as  the  viscous  cutoff  region. 

At  intermediate  wavenumbers,  the  exponential  growth  term  dominates  and  the 
equilibrium  spectrum  predicted  by  this  model  is  larger  than  the  Pierson-Moskowitz 
spectrum.  The  amount  of  this  enhancement  depends  on  the  wind  speed  and  also  on 
the  value  of  the  constant  y .  An  example  plot  of  the  equilibrium  wave  height 
spectrum  for  y=1  ,  using  the  growth  rate  suggested  by  Plant  (1982),  is  shown  in 
Figure  1.  The  wavenumber  dependence  shown  here  is  qualitatively  similar  to  that 
suggested  by  Pierson  and  Stacy  (1973)  and  Bjerkaas  and  Riedel  (1979).  However, 
the  detailed  structure  shown  in  the  measurements  by  Jahne  and  Riemer  (1990)  is  not 
accounted  for  by  this  model  (see  sections  2.5  and  4.2  below).  Multiplying  the  height 
spectrum  by  Jc2  cos2  {<)>-<}>„)  and  k2  sin2  (4>-4>„)  and  integrating  over  all 
wavenumbers  yields  the  slope  variances  in  the  upwind  and  crosswind  directions, 
respectively,  as  shown  in  Figure  2.  These  slope  variances  agree  fairly  well  with 
those  measured  by  Cox  and  Munk  (1954)  which  are  also  shown  in  Figure  2. 
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LEGEND 

Up-wind  component _ 

Cross-wind  component 
a  Cox  Sc  Munk  dato  (upwind) 

+  Cox  Sc  Munk  dato  (crosswind) 

Figure  2.  Comparison  of  Slope  Variances  Calculated  From  Spectral  Balance 
Model  With  Measurements  by  Cox  and  Munk  (1954). 
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2.2.  WAVE-CURRENT  INTERACTIONS 

The  net  source  function  discussed  in  the  previous  section  can  be  rewritten  in 
terms  of  the  action  spectral  density  JV(k,4>)  =p cS(k,4>)  as 

-^=oc/  +  ($-4vk2)  N-y'(i}k*N2 
at 

where  a/=pca  and  y'=y /pc.  This  equation  can  be  used,  along  with  the  ray 
equations 


+u 

-jr  —  -k  & 

dt  ** 

dt 

Xdx  ydx 

dy=  c  +v 

dky_ 

dt  ** 

dt 

xdy  y  By 

to  calculate  the  evolution  of  a  wave  packet  as  it  interacts  with  a  variable  surface 
current.  Alternatively,  these  five  equations  can  be  combined  into  a  single  partial 
differential  equation  as  discussed  in  Lyzenga  and  Bennett  (1988).  The  paper 
reproduced  in  Appendix  A  discusses  an  analytic  solution  of  this  equation,  obtained  by 
approximating  the  net  source  function  by  -pr  (N-N0)  where  Na  is  the  equilibrium 
action  spectrum  and  Pr  is  a  relaxation  rate.  Using  a  Taylor  expansion  of  the  net 
source  function  described  above,  this  relaxation  rate  may  be  written  as 

Pr=v^(P~4vic2)i-*-4aY<i)ic4  . 

Note  that  Pr~  P  in  the  intermediate  wavenumber  range  where  the  exponential  growth 
term  dominates,  and  Pr=4vJc2  in  the  viscous  cutoff  region. 

An  example  calculation  of  the  change  in  the  wave  spectrum  due  to  a 
converging  current  is  shown  in  Figure  3.  Additional  examples  are  discussed  in 
Appendix  A  and  in  section  4.2  below.  A  general  characteristic  of  these  calculations 
is  that  the  maximum  effects  are  predicted  at  wavelengths  on  the  order  of  1  meter,  and 
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1  m 


"10  cm' 


-J  =  —0.01  sec 


\dx) 


=  0 


100  m 


U  =  5  m/s 


Contour  Plot  of  the  Change  in  the  Curvature  Spectrum,  i.e. 
k*S( k,  4>)  -k*Sa  (k,  <J>)  ,  for  a  Converging  Current.  Fractional 
Change  in  the  Spectral  Density  is  on  the  Order  of  Unity  at  X=1  m  and 
One  Percent  at  X=1  cm  (Note  Logarithmic  Wavenumber  Scale). 
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very  small  changes  in  the  wave  spectrum  are  predicted  at  wavelengths  of  less  than  10 
cm.  These  predictions  are  apparently  at  odds  with  the  observation  of  large  changes  in 
microwave  backscatter  associated  with  current  gradients  on  the  ocean  surface,  since 
the  microwave  backscatter  is  generally  considered  to  be  governed  by  the  surface  wave 
spectral  density  at  wavelengths  comparable  to  the  radar  wavelength.  The  studies 
discussed  in  the  remainder  of  this  report  were  motivated  in  large  part  by  this  apparent 
discrepancy.  A  more  quantitative  evaluation  of  the  accuracy  of  the  predictive  models 
will  be  presented  in  section  4  below. 


2.3.  SURFACTANT  EFFECTS 

The  viscous  damping  of  surface  waves  is  enhanced  by  the  presence  of 
surfactant  films.  A  linearized  theory  of  this  damping  effect  was  developed  by 
Dorrestein  (1951)  and  Levich  (1962).  Empirical  validation  of  this  theory  was 
provided  by  Cini  and  Lombardini  (1981).  The  main  parameter  influencing  the 
amount  of  wave  damping  is  the  surface  elasticity  E  which  describes  the  change  in 
surface  tension  upon  compression  of  the  film.  Diffusion  or  adsorption  effects  can  be 
incorporated  by  considering  this  to  be  a  complex  quantity.  However,  in  the  simpler 
case  in  which  the  elasticity  is  real,  the  surface  wave  damping  rate  is  given  by 


0d=4vJc2 


l-X+XY 

1-2X+2X2 


where 


X-  g*2- 
PWv/2vq 


and 


Y _  Ek 
4pvo) 


(Cini  et  al. ,  1987).  Figure  4  shows  a  plot  of  the  viscous  damping  rate  for  a  clean 
surface  (E=  0)  and  for  a  surface  covered  with  oleyl  alcohol,  which  has  an  elasticity 
E= 22.5  dyne/cm.  The  damping  rate  for  the  film-covered  surface  has  a  characteristic 
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peak  at  a  wavenumber  of  about  2  rad/cm,  which  is  caused  by  a  resonant  coupling  of 
the  surface  wave  energy  into  longitudinal  or  Marangoni  waves  within  the  surface  film. 

This  damping  rate  can  be  used  in  the  net  source  function,  as  discussed  in 
Section  2.1,  in  order  to  determine  the  change  in  the  equilibrium  spectrum  due  to  the 
presence  of  a  surface  film.  Example  results  of  such  a  calculation  are  shown  in  Figure 
5  for  an  oleyl  alcohol  film. 

2.4.  ATMOSPHERIC  STABILITY  EFFECTS 

The  wind- wave  growth  rate  introduced  by  Plant  (1982)  involves  the  friction 
velocity  u.=/c7p  where  x  is  the  wind  stress  and  p  is  the  density  of  air.  The  wind 
stress  is  dependent  on  atmospheric  stability  conditions  as  well  as  the  wind  speed.  An 
empirical  relationship  among  these  variables,  as  determined  from  a  set  of  214  records 
over  a  range  of  windspeeds  from  2  to  21  m/sec  (Geemaert,  1990),  can  be  written  as 

uf=  (2.58  U+  0 . 4  9  £/2  +  0 . 07  £73  - 1 . 06  A T)  x  10'3  m2/sec2 

where  U  is  the  windspeed  (in  m/sec)  measured  at  10  m  above  the  surface  and  A  T  is 
the  difference  between  the  air  and  water  temperatures  (in  °  c).  Using  this  relationship 
to  calculate  the  growth  rate,  and  solving  for  the  equilibrium  spectrum  as  discussed  in 
section  2.1,  yields  the  results  shown  in  Figure  6  for  [7=3  m/sec  and  two  values  of 
A  T  corresponding  to  stable  ( A  T  >  0 )  and  unstable  ( A  T  <  0 )  conditions.  Note  that  the 
largest  effects  are  predicted  at  very  short  wavelengths. 

2.5.  WAVE  BREAKING  EFFECTS 

The  net  source  function  discussed  in  section  2. 1  includes  a  dissipation  term 
which  is  thought  to  be  due  primarily  to  the  effects  of  wave  breaking  (Phillips,  1977, 

1 985).  However,  it  does  not  incorporate  any  corresponding  source  of  wave  energy  or 
spectral  density.  There  is  ample  empirical  evidence  that  small-scale  roughness  is 
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U=  3.0  m/s,  Td=  0.0  C,  E=  0.0  mN/m 
U=  3.0  m/s,  Td=  0.0  C,  E=  22.5_mN/m 


Figure  5.  Equilibrium  Wave  Spectrum  Calculated  From  Spectral  Balance  Model 
for  Clean  Water  and  for  a  Surface  Covered  With  Oleyl  Alcohol 
Monolayer. 
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generated  by  breaking  waves  (e.g.  Banner  and  Fooks,  1985),  although  there  is 
consider.  ’:!e  uncertainty  as  to  the  details  of  this  process.  This  section  describes  a  first 
attempt  to  model  the  effects  of  wave  breaking  on  the  surface  roughness.  The  effects 
of  this  change  in  roughness  on  the  radar  backscatter  are  considered  in  section  3.3 
below. 

A  procedure  for  predicting  the  fraction  of  the  surface  covered  by  breaking 
wave  crests  has  been  developed  and  to  some  extent  validated  by  several  authors.  The 
procedure  assumes  that  there  is  a  threshold  on  the  downward  vertical  acceleration  at 
the  surface  of  the  water,  beyond  which  the  surface  becomes  unstable  and  breaks.  A 
breaking  threshold  of  -g! 2  was  suggested  by  Longuet-Higgins  (1969),  this  being  the 
limiting  acceleration  near  the  crest  of  a  Stokes  wave.  Laboratory  observations  of 
wave  breaking  by  Ochi  and  Tsai  (1983)  were  shown  by  Srokosz  (1983)  to  be 
consistent  with  a  threshold  of  -ag,  where  a  *0.4.  Using  this  threshold  and  assuming 
that  the  values  of  the  vertical  acceleration  are  normally  distributed  with  mean  zero 
and  standard  deviation  oa,  as  in  Snyder  and  Kennedy  (1983),  the  fraction  of  the 
surface  covered  by  breaking  water  is  given  by 
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exp 
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For  gravity  waves,  the  variance  of  the  vertical  acceleration  is  given  by 
wc  2n 

o\  =  j  v>*Su  ({*>)  do>  =  J  J  g2k2S(k,$)  kdkdfy  =  g2o2g 
0  0  0 

where  o|  is  the  slope  variance  and  kc  is  a  cutoff  wavenumber  which  defines  the 
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spatial  scale  of  the  breaking  waves.  A  plot  of  the  breaking  fraction  as  a  function  of 
the  slope  variance  is  shown  in  Figure  7.  The  cutoff  wavenumber  is  somewhat 
arbitrary,  but  for  the  present  purposes  it  is  assumed  to  lie  near  the  boundary  between 
gravity  and  capillary  waves.  The  slope  variance  in  this  equation  may  therefore  be 
associated  with  that  measured  by  Cox  and  Munk  (1954)  in  the  presence  of  a  surface 
slick,  which  would  presumably  damp  out  the  capillary  waves  as  discussed  in  section 
2.3  above.  The  total  (upwind  plus  crosswind)  slope  variance  under  these  conditions 
was  found  to  be  related  to  the  wind  speed  by  the  equation 

o|  =  ol  +  02u  =  0. 008  +  0. 00156W 

where  W  is  the  wind  speed  (in  m/sec)  recorded  at  41  feet  above  sea  level.  The 
fractional  breaking  area  using  this  equation  for  the  slope  variance  is  plotted  in  Figure 
8.  For  wind  speeds  between  3  and  15  m/s  the  breaking  fraction  is  approximated  quite 
closely  by  the  equation 

fb~Q  .  0044|-^|2  57  3<W<15m/s. 

This  is  roughly  half  of  the  fractional  whitecap  coverage  as  given  by  Spillane  et  al 
(1986),  although  is  is  somewhat  closer  to  the  observations  of  Kondo  et  al  (1973). 

Given  the  foregoing  macroscopic  description  of  the  surface  in  terms  of  the 
breaking  fraction  fb,  the  next  step  is  to  attempt  to  statistically  describe  the 
microstructure  within  the  breaking  regions.  Here  we  are  faced  with  the  problem  that 
wave  breaking  may  encompass  a  wide  variety  of  phenomena,  ranging  from  parasitic 
capillary  generation  at  the  onset  of  small-scale  breaking  to  the  formation  of  droplets, 
bubbles,  and  foam  during  the  most  energetic  large-scale  breaking  events.  Although 
the  latter  cases  are  the  most  obvious  and  easily  observed,  the  breaking  criterion 
discussed  above  suggests  that  a  large  fraction  of  the  breaking  surface  area  may  be 
associated  with  relatively  short  waves.  For  a  k'4  spectrum,  the  slope  variance  within 
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the  wavenumber  interval  from  k1  to  k2  is  proportional  to  In  {k2/k1)  .  Thus,  for 
example,  wavelengths  from  10  cm  to  1  meter  contribute  the  same  amount  as 
wavelengths  from  1  to  10  meters.  Therefore,  features  associated  with  small-scale 
breaking  waves,  such  as  parasitic  capillary  wave  generation,  may  be  as  important  as 
large-scale  breaking  events  for  the  purposes  considered  here. 

Laboratory  measurements  by  Banner  and  Fooks  (1985)  indicate  that  the 
roughness  generated  by  20-30  cm  breaking  waves  is  fairly  narrow-banded.  Field 
measurements  by  Kondo  et  al  (1973)  associated  breaking  waves  with  local  increases 
in  the  height  variance  calculated  over  time  intervals  of  0.2  sec.  Typical  values  of  the 
r.m.s.  surface  height  for  these  high-frequency  components  were  on  the  order  of  0.5-1 
cm  in  the  presence  of  breaking  waves.  Values  of  this  r.m.s.  height  were  also 
calculated  for  all  observations  (with  and  without  breaking)  and  plotted  versus  the  wind 
speed. 

For  a  given  measurement  interval,  the  height  variance  associated  with  the 
roughness  caused  by  wave  breaking  may  be  assumed  to  ’  e  proportional  to  the  fraction 
of  the  surface  covered  by  breaking  water  during  this  interval,  and  may  be  written  as 

hl=£bh2b 

where  hi  is  the  height  variance  within  the  breaking  regions.  The  spatial  structure  of 
both  breaking  and  non-breaking  regions  can  be  described  by  means  of  the 
autocovariance  function 

<j>(r)  =(ii  (rOt)  (r'+r)) 

where  q  (r)  is  the  surface  height  and  the  brackets  indicate  ensemble  averages.  It 
will  be  assumed  here  that  the  surface  within  a  breaking  region  may  be  characterized 
by  a  covariance  function  having  a  correlation  length  zc  which  is  much  smaller  than 
the  correlation  length  in  the  absence  of  breaking. 
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Assuming  the  small-scale  height  variations  within  the  breaking  regions  are 
uncorrelated  with  the  larger-scale  height  variations  on  the  underlying  surface,  the  total 
height  autocovarianca-  function  can  be  written  as 

4>(r)  =  <Mr)  +  ^(r)  . 

It  is  further  assumed  that  the  shape  of  the  breaking  wave  contribution  to  the 
covariance  function  is  approximately  gaussian  and  isotropic,  and  can  thus  be  written 
as 

4>b(r)  =  fbhl  e~zl/z* 

where  rc  is  the  correlation  length  within  the  breaking  region. 

Using  this  expression  for  the  surface  height  autocovariance  function,  the 
corresponding  height  spectrum  can  be  written  as 

S(k.Q)  =  - — —  f  fi>  (r)  ei*rco8  (fl-fl')  z  dr  jQf 

(2n)ZJ  J 

-  Jo(kz)idz  *  s^ik.e) 

Detailed  measurements  of  the  short  wavelength  portion  of  the  wave  spectrum 
in  a  large  wind  wave  facility  have  been  made  by  Jahne  and  Riemer  (1990).  These 
measurements  show  a  weak  dependence  of  the  spectrum  on  the  wind  speed  or  friction 
velocity  at  wavenumbers  below  1  rad/cm  and  a  much  stronger  dependence  at  higher 
wavenumbers,  with  a  sharp  cutoff  in  the  spectrum  at  a  wavenumber  of  approximately 
8  rad/cm.  The  wind  speed  dependence  for  wavenumbers  between  1  and  8  rad/cm  and 
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friction  velocities  between  about  10  and  50  cm/s  is  described  as  approximately  a 
power  law  in  the  friction  velocity,  with  an  exponent  between  2.5  and  3.  For  the 
lowest  wind  speeds  the  dimensionless  curvature  spectrum  or  "degree  of  saturation” 
B(k),  which  is  defined  as  If  times  the  height  spectrum,  falls  off  rapidly  between 
wavenumbers  of  1  and  2  rad/cm. 

One  possible  interpretation  of  Jahne  and  Riemer’s  measurements  is  that  for 
wavenumbers  below  1  rad/cm,  the  spectrum  results  from  a  balance  between  the  wind 
input  and  various  dissipation  mechanisms  as  discussed  in  section  2.1  above,  whereas 
for  ^ner  wavenumbers  the  spectrum  is  dominated  by  wave  breaking  effects.  Under 
this  hypothesis,  the  rapid  increase  in  the  spectrum  with  wind  speed  at  these 
wavenumbers  is  then  due  to  the  strong  dependence  of  the  breaking  fraction  fb  on  the 
slope  variance.  The  observed  cutoff  in  the  spectrum  at  8  rad/cm  is  consistent  with  a 
value  of  re  of  roughly  0.5  cm.  Furthermore,  the  observed  spectral  levels  in  the 
region  of  1  to  8  rad/cm  appear  to  be  predicted  fairly  well  by  assuming  a  value  of 
hb~rc  and  using  the  breaking  fraction  calculated  as  described  earlier  in  this  section. 
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3.0  RADAR  BACKSCATTER  MODELING 

The  problem  of  predicting  the  statistics  of  the  radar  backscatter  from  the  ocean 
surface  is  one  of  long  standing.  There  are  well-known  solutions  of  this  problem  in 
the  limiting  cases  where  the  electromagnetic  wavelength  is  either  much  longer  or 
much  shorter  than  the  relevant  spatial  scales  of  the  surface  roughness.  However,  for 
electromagnetic  radiation  in  the  microwave  region  of  the  spectrum,  and  at 
intermediate  angles  of  incidence,  the  spatial  scales  for  typical  ocean  surfaces  are  such 
that  neither  of  these  limiting  solutions  is  valid. 

The  most  widely  accepted  approximate  solution  within  this  intermediate 
scattering  regime  is  the  composite  or  two-scale  model.  Some  problems  associated 
with  this  model  are  discussed  in  section  3.2  below.  An  alternative  approach  is 
provided  by  the  Kirchhoff  or  physical  optics  approximation.  This  approximation  is 
reviewed  in  the  following  section,  and  is  applied  to  the  problem  of  breaking  wave 
backscatter  in  section  3.3. 


3.1.  KIRCHHOFF  APPROXIMATION 

Electromagnetic  radiation  striking  the  surface  of  a  highly  conducting  medium, 
such  as  the  ocean,  induces  a  surface  current  Jm  (r)  .  Once  this  current  is  known,  the 
scattered  magnetic  field  can  be  calculated  from  the  equation 

Ht (r)  =JVG{jr-r')  xJB(r')  dx' 

where 


G(r-r/) 


e  lk\x-z'\ 

4n  |r-r'| 


The  Kirchhoff  or  physical  optics  approximation  is  obtained  by  assuming  that  the 
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surface  current  is  given  by 


J,(r)  =2 air)  xE1lz) 

where  n(r)  is  the  surface  normal  and  Bi  (r)  is  the  incident  magnetic  field.  This  is 
equivalent  to  the  current  which  would  be  present  on  an  infinite  flat  plane  tangent  to 
the  surface  at  the  location  r.  If  the  surface  height  r\  {x,y)  is  assumed  to  be  a 
Gaussian-distributed  random  variable,  the  radar  cross  section  per  unit  area  is  then 
given  by 

a0=-~klsec2Qj jr(x,y)  e~2ikb* dxdy 

where  kQ  is  the  radar  wavenumber,  6  is  the  incidence  angle,  kh~ k0sinQ , 
kz=kacosQ  and 

T(x,y)  =exp{4Jtf  [^(x,y)  -6(0,0)  ] } 

where  6  (x,y)  is  the  surface  height  autocovariance  function  (Holliday  et  al.,  1986). 

In  the  limit  as  kz~<*>,  T(x,y)  is  non-zero  only  in  an  infinitesimal  region 
around  the  origin,  over  which  the  covariance  function  6  (*,  y)  can  be  expanded  in  a 
Taylor  series  to  second  order  in  x  and  y.  The  integral  for  o0  can  then  be  evaluated 
analytically  to  yield  the  well-known  specular  scattering  result  (Barrick,  1968).  On  the 
other  hand,  as  kz-*  0  the  small-argument  approximation  can  be  used  for  the 
exponential  function  in  the  above  equation,  with  the  result  that  o a  becomes 
proportional  to  the  surface  height  spectral  density  at  the  "Bragg"  wavenumber 
kB=2kh=2k0sind.  However,  the  constant  of  proportionality  in  the  resulting 
expression  is  different  from  that  obtained  by  the  small  perturbation  method  (SPM).  It 
has  been  shown  by  Holliday  (1987)  that  if  a  second  iteration  is  performed  by 
calculating  the  current  induced  by  the  first-order  scattered  field  in  the  Kirchhoff 
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approximation,  the  SPM  result  is  reproduced  for  both  vertical  and  horizontal 
polarization. 

It  can  be  seen,  therefore,  that  although  the  Kirchhoff  approximation  is  inferior 
to  the  SPM  in  some  respects,  it  has  certain  advantages  in  that  it  does  not  rely  on  the 
assumption  that  the  waveheight  is  much  smaller  than  the  electromagnetic  wavelength. 
This  makes  the  approximation  attractive  for  some  cases,  such  as  the  one  discussed  in 
section  3.3,  where  the  wave  height  may  not  be  small  enough  for  the  SPM  to  be  valid 
and  for  which  the  two-scale  model  is  not  applicable  because  the  roughness  is  basically 
on  a  single  length  scale. 

3.2.  TWO-SCALE  MODEL 

There  are  several  variations  and  several  methods  of  deriving  the  two-scale 
model,  but  the  basic  idea  is  to  apply  the  small  perturbation  or  Bragg  scattering  model 
to  small  patches  of  the  surface,  each  of  these  patches  having  a  different  slope  or 
orientation  with  respect  to  the  overall  scattering  surface.  The  total  scattered  power  is 
then  obtained  by  summing  or  integrating  over  all  possible  slopes  of  the  Bragg 
scattering  patches.  There  is  some  ambiguity  regarding  the  scale  separation 
wavelength,  which  influences  the  distribution  of  slopes  included  in  the  integration, 
and  also  regarding  the  method  of  performing  this  integration.  In  addition,  the  two- 
scale  model  makes  certain  assumptions  about  the  surface  height  and  slope  statistics 
(namely,  that  they  are  uncorrelated)  which  may  be  violated  at  small  scales  and  in  the 
presence  of  wave  breaking  and  parasitic  capillary  wave  generation. 

The  radar  cross  section  given  by  the  two-scale  model  can  be  written  in  its 
simplest  form  as 

oa=ffoB(r)xti)y)p(rtxrr)y)  dr)xdr)y 

where  ob(tix,  rjy)  represents  the  Bragg  scattering  cross  section  per  unit  area  for  a 
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surface  element  having  slope  components  tix  and  r\y,  and  p  {t}x,  tjy)  represents  the 
probability  density  function  for  these  slopes.  This  integral  can  be  evaluated 
numerically,  using  well-known  expressions  for  the  Bragg  scattering  cross  section  as 
described  for  example  in  Lyzenga  and  Bennett  (1988).  An  example  of  the  results  for 
vertically  and  horizontally  polarized  radiation,  assuming  a  k*  short  wave  spectrum 
and  a  45°  incidence  angle,  are  shown  in  Figures  9  and  10  as  a  function  of  the  slope 
variance  <r,2. 

Alternatively,  the  Bragg  scattering  cross  section  can  be  expanded  as  a  Taylor 
series  in  the  slopes  rix  and  i\y  and  the  integration  can  be  carried  out  analytically.  In 
fact,  it  is  argued  by  Plant  (1986)  that  it  is  "inconsistent"  to  carry  out  this  expansion 
beyond  second  order  in  the  surface  slope.  This  implies  that  the  radar  cross  section  is 
a  linear  function  of  the  slope  variance,  as  shown  by  the  dashed  lines  in  Figures  9  and 
10.  Obviously,  there  is  a  large  difference  between  the  results  for  horizontal 
polarization  when  the  slope  variance  exceeds  about  0.03,  or  the  r.m.s.  tilt  angle 
exceeds  10°. 

3.3.  BACKSCATTER  FROM  BREAKING  WAVES 

Measurements  of  the  X-band  radar  backscatter  from  breaking  waves  at  near¬ 
grazing  angles  of  incidence  led  Lewis  and  Olin  (1980)  to  propose  that  breaking 
regions  may  be  considered  as  perfect  isotropic  reflectors,  with  all  of  the  incident 
energy  being  scattered  uniformly  throughout  the  upper  hemisphere.  This  would  imply 
a  radar  cross  section  per  unit  area  on  the  order  of  unity,  or  more  precisely 

o0~2 cos© 

where  0  is  the  angle  of  incidence.  Field  measurements  of  microwave  backscatter 
made  at  intermediate  incidence  angles  have  also  been  interpreted  in  terms  of  the 
contributions  from  wave  breaking  by  Jessup  et  al  (1990,  1991a,  1991b)  as  discussed 
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RADAR  CROSS  SECTION 


VERTICAL  POLARIZATION 


Figure  9.  Radar  Cross  Section  Versus  Large-Scale  Slope  Variance  Calculated 

From  Two  Versions  of  the  Two-Scale  Model,  for  Vertical  Polarization 
at  45°  Incidence.  Solid  Line  Indicates  Result  of  a  Numerical 
Integration  Over  Wave  Slopes,  as  Described  by  Lyzenga  and  Bennett 
(1988),  and  Dashed  Line  Indicates  Analytical  Result  Obtained  by  Plant 
(1986). 
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HORIZONTAL  POLARIZATION 


SLOPE  VARIANCE 


Figure  10.  Radar  Cross  Section  Versus  Large-Scale  Slope  Variance  Calculated 
From  Two  Versions  of  the  Two-Scale  Model,  for  Horizontal 
Polarization  at  45°  Incidence.  Solid  Line  Indicates  Result  of  a 
Numerical  Integration  Over  Wave  Slopes,  as  Described  by  Lyzenga 
and  Bennett  (1988),  and  Dashed  Line  Indicates  Analytical  Result 
Obtained  by  Plant  (1986). 


j>ERIM _ 

below.  Additional  laboratory  measurements  of  the  radar  backscatter  from  breaking 
waves  have  been  made  by  Banner  and  Fooks  (1985),  Melville  et  al.  (1988),  and 
Trizna  et  al.  (1991). 

Several  attempts  have  been  made  to  model  the  backscatter  from  breaking 
waves  deterministically,  based  on  various  assumptions  about  the  geometric  structure 
of  the  surface.  Wetzel  (1981,  1990)  has  constructed  models  for  breaking  wave 
backscatter  based  on  assumptions  about  the  detailed  shapes  of  the  scattering  elements 
within  an  actively  breaking  region.  Lyzenga  et  al.{ 1983)  proposed  a  wedge 
diffraction  model  to  account  for  scattering  from  the  sharp  crests  of  breaking  waves, 
and  Kwoh  and  Lake  (1984)  made  laboratory  measurements  and  calculations  to 
evaluate  the  relative  effects  of  the  sharp  crest  and  the  parasitic  capillary  waves 
associated  with  a  finite-amplitude  surface  gravity  wave. 

Most  statistical  scattering  models  of  the  microwave  backscatter  from  the  ocean 
surface  at  intermediate  incidence  angles  have  been  based  on  Bragg  or  two-scale 
approximations  (e.g.  Valenzuela,  1978).  Phillips  (1988)  has  suggested  that  the  radar 
cross  section  of  the  ocean  surface  may  be  represented  as  the  sum  of  separate 
contributions  from  Bragg  scattering  and  individual  wave  breaking  events  but  did  not 
specify  the  scattering  mechanism  for  these  events. 

In  the  following,  the  Kirchhoff  approximation  discussed  in  section  3.1  is 
combined  with  the  statistical  description  of  a  breaking  wave  surface  discussed  in 
section  2.5  in  order  to  construct  a  model  for  the  backscatter  from  breaking  waves. 

Using  the  Gaussian  surface  height  autocovariance  function  discussed  in  section 
2.5,  the  function  T(x,y)  which  enters  the  Kirchhoff  approximation  can  be  written  as 

ru)=rb(r)rnbu) 

=  crJb(r)  -  rto(oo) ]  r^(zr)  *■  iy^r^r) 

where 
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T*(r)  =  exp{4*|C4»J>(r)  -<M°>  ]},  Tb(<») 

and 

^(r)  =  exp{4^  [4>^(r) -^(0)  3 }  . 

If  we  further  assume  that  the  correlation  length  in  the  absence  of  breaking  is  much 
longer  than  rc,  we  can  approximate  this  function  as 

r(r)*[r2>{r)-ri(~)]rnA(o)  + 

-iyri-iyo.)  ♦  r^oojr^tr) . 


Substituting  this  into  the  Kirchhoff  approximation,  the  radar  cross  section  per  unit 
area  in  the  presence  of  wave  breaking  becomes 


<V=  ab+°i>be 


-Aklhl 


where  is  the  radar  cross  section  per  unit  area  in  the  absence  of  breaking,  and 
ob= A  kl  sec2e|J  [Tb  ( r )  -rb  ( »)  ]  e  ~2ik»rcoa*'r  dr  d4,/ . 


To  evaluate  this  expression,  we  make  another  simplifying  assumption  and  approximate 
the  integrand  by  a  second  gaussian  function,  i.e. 

Tb{r)  -iy«)  «(l-e'4*’**)e~ri/£j 

where  L  is  equal  to  re  in  the  low-frequency  limit  and  approaches  the  value 

ic!  (2kzhb)  in  the  high-frequency  limit.  In  the  intermediate  frequency  region,  the 

value  of  L  may  be  approximated  by  the  formula 
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Am—m* 

Hie  Kirchhoff  integral  can  then  be  readily  evaluated  to  yield  the  result 

ab=kl  L2 e'k*L* (l-e~ik‘b*) sec2Q 


Before  comparing  this  expression  with  observations,  some  of  its  properties 
may  be  noted.  First,  it  can  be  seen  that  in  the  high-frequency  limit  the  radar  cross 
section  approaches  the  value 


ob= 


as  kQ—  <*> 


which  is  recognized  as  the  familiar  specular  scattering  result.  Second,  in  the  low- 
frequency  limit  or  for  small  grazing  angles,  the  small-argument  approximation  can  be 
used  for  the  exponential  function,  and  the  cross  section  approaches  the  limiting  value 


oJb=4iCoh|r|e'^r*  as  k0~  0  or  0-it/2 . 

This  is  equivalent  to  the  result  obtained  by  using  the  small  perturbation  method 
(SPM),  except  that  in  this  case  ob  is  multiplied  by  cos40  for  horizontal  polarization 
and  (l+sin40)  for  vertical  polarization.  Plots  are  shown  in  Figure  11  of  the 
backscatter  cross  section  calculated  from  the  SPM,  from  a  numerical  integration  of 
the  Kirchhoff  integral,  and  from  the  analytical  approximation  to  the  Kirchhoff  integral 
developed  above.  Although  the  Kirchhoff  model  does  not  predict  the  polarizati  m 
dependence  of  the  backscatter,  it  does  at  least  provide  a  transition  from  the  low- 
frequency  to  the  high-frequency  cases  and  may  therefore  be  expected  to  be  more 
accurate  than  the  SPM  in  the  intermediate  frequency  case,  although  the  limits  of 
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Figure  11. 
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Comparison  of  X-Band  Radar  Backscatter  Cross  Sections  Calculated 
From  Kirchhoff  and  SPM  Models  for  Surface  Described  by  a 
Gaussian  Correlation  Function  With  rms  Height  and  Correlation 
Length  Both  Equal  to  0.5  cm. 
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validity  of  both  of  these  models  are  still  uncertain. 

Using  the  approximate  Kirchhoff  expression  for  ab  derived  above,  and 
assuming  rc=hb=0.5  cm  as  discussed  in  section  4,  the  radar  cio^s  section  per  unit 
area  for  k0—2  rad/cm  (X-band)  is  equal  to  approximately  1.5  at  grazing  incidence, 
which  agrees  with  observations  by  Lewis  and  Olin  (1980).  For  the  case  of  large-area 
averages  where  a  fraction  of  the  surface  is  covered  by  breaking  water,  the 
contribution  of  breaking  waves  to  the  backscatter  can  be  written  as 

where fb  is  the  breaking  fraction  discussed  in  section  2.5.  This  expression  may  be 
compared  with  the  observations  of  Jessup  et  a/.  (1991a,  1991b)  by  using  the  breaking 
fraction  fb  as  a  function  of  wind  speed  as  shown  in  Figure  8,  and  the  values  of  rc  and 
hb  as  estimated  in  section  2.5.  For  a  radar  wavenumber  of  2.94  rad/cm  (K„  band)  and 
an  incidence  angle  of  45°,  the  radar  cross  section  per  unit  area  for  the  breaking  areas 
is  approximately  0.66  and  the  average  contribution  of  these  areas  to  the  entire  data  set 
is  equal  to  0.66  fb.  This  contribution  is  plotted  versus  the  friction  velocity  (assumed 
to  be  1/30  of  the  wind  speed)  and  compared  with  measurements  from  Jessup  et 
al.  (1991b)  in  Figure  12. 
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Figure  12.  Comparison  of  Model  Predictions  (Dashed  Line)  With  Measurements 
by  Jessup  et  al  (1991b)  of  Breaking  Wave  Contributions  to  the 
Average  Radar  Cross  Section  of  the  Ocean  Surface. 
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4.0  MODEL  VALIDATION  STUDIES 

Aside  from  comparisons  with  previously  published  data,  model  validation 
activities  undertaken  during  this  project  centered  on  the  ONR/NRL  High  Resolution 
field  experiment.  This  activity  involved  three  areas  of  participation:  (1)  experiment 
planning  and  pre-experiment  model  predictions,  (2)  deployment  of  a  Doppler  radar  on 
the  Woods  Hole  LADAS  platform,  and  (3)  analysis  of  SAR  and  in  situ  data,  and 
comparison  with  model  predictions. 

A  number  of  interesting  features  associated  with  the  edge  of  the  Gulf  Stream 
were  observed  in  both  the  SAR  and  RAR  images  collected  during  the  first  High 
Resolution  experiment.  These  include  (1)  numerous  bright  lines,  some  aligned 
parallel  to  the  Gulf  Stream  boundary  and  some  oriented  at  large  angles  to  this 
boundary,  (2)  dark  narrow  lines  parallel  to  the  Gulf  Stream  currents,  and  (3)  an 
overall  higher  return  over  the  Gulf  Stream  as  compared  with  the  shelf  water  to  the 
west  of  the  stream. 

The  features  in  the  first  category  are  presumably  caused  by  strong  shearing 
and  converging  currents.  Although  these  currents  were  not  always  clearly  observed 
in  the  ADCP  data  because  of  their  shallow  depth  and  rapid  spatial  variations,  some 
information  has  been  obtained  by  analyzing  a  sequence  of  RAR  images  (F.  Askari, 
personal  communication).  An  analysis  of  the  radar  backscatter  variations  associated 
with  these  features  is  presented  in  section  4.2  below. 

The  second  category  of  features  may  be  due  to  the  accumulation  of  surfactant 
materials  along  weaker  convergence  or  shear  regions.  Evidence  for  the  association  of 
such  features  with  shearing  currents  has  been  observed  during  previous  experiments  in 
which  ship  wakes  crossing  similar  dark  lines  were  observed  to  be  displaced  on  either 
side  of  the  lines  (A.  Ochadlick,  personal  communication).  However,  the  connection 
between  the  dark  lines  observed  along  the  edge  of  the  Gulf  Stream  and  and  any  sharp 
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current  gradients  in  this  region  has  not  yet  been  established. 

The  higher  radar  returns  observed  over  the  Gulf  Stream  in  some  of  the  radar 
images  may  be  caused  by  atmospheric  stability  effects  associated  with  the  higher 
water  temperatures  in  the  Gulf  Stream,  as  discussed  in  section  2.4  of  this  report. 
Further  testing  of  this  hypothesis  is  required  using  the  meteorological  data  collected 
during  the  pilot  experiment. 

Several  of  the  features  observed  in  the  airborne  radar  images  collected  during 
the  pilot  experiment  have  been  identified  in  the  Doppler  radar  data  set  as  well.  The 
intention  was  to  compare  this  data  with  in  situ  measurements  of  the  surface 
roughness,  surface  currents  and  wind  stress  in  order  to  test  the  wave  spectrum  and 
radar  backscatter  models  described  above.  Activities  carried  out  in  pursuit  of  this 
goal  are  described  in  the  following  section. 

4.1.  DOPPLER  RADAR  MEASUREMENTS 

During  this  project,  a  simple  Doppler  radar  system  was  assembled  using  a 
commercial  motion  sensing  device  built  by  Microwave  Associates,  and  this  system 
was  deployed  on  the  Woods  Hole  LADAS  catamaran  during  the  High  Resolution  field 
experiment.  The  MA86735  device  consists  of  a  Gunn  diode  source  in  a  tuned  cavity, 
which  is  coupled  to  a  short  section  of  waveguide  in  which  two  Schottky  diode  mixers 
are  also  mounted  (see  Figure  13).  The  other  end  of  the  waveguide  is  connected  to  a 
horn  antenna  which  radiates  the  microwave  signal  and  receives  the  radiation 
backscattered  by  objects  or  surfaces  within  its  field  of  view.  The  electric  field  at  each 
of  the  detectors  can  be  represented  by 

Ei  =  Eti  +  Eii 

where  Eti  is  the  outgoing  or  transmitted  field  and  Ezl  is  the  incoming  or  received 
field  at  detector  i.  If  this  detector  is  located  a  distance  x,  from  the  source,  the 
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Frequency  =  10.525  GHz 
Polarization  =  Vertical 
Antenna  Beamwidth  =  12  deg. 


Figure  13.  Schematic  Diagram  of  Doppler  Radar.  Gunn  Diode  Source  is  on  the 
Left.  Vdl  and  Vc  are  the  In-Phase  and  Quadrature  Signals  From  the 
Detector  Diodes. 
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transmitted  field  is  given  by 


Eti 


and  the  received  field,  for  a  scatterer  located  at  a  distance  x  from  the  source  is 


-ti 


g  lej(2kx-kxi-at) 


where  k-2-n/k  is  the  radar  wavenumber,  u=2tc/  is  the  radar  frequency  in 
radians/sec  and  r  is  the  ratio  of  the  transmitted  to  received  field  strength  (which 
depends  on  the  radar  cross  section  of  the  scatterer,  its  range  and  position  within  the 
antenna  beam).  Assuming  that  r<l,  the  voltage  on  the  detector  is  then  given  by 

Vi  =  | Eci+Eri\  =  E0\l+rej2kl*~Xl)  |  »  B0ll+zcoe(2kx-+i)) 

where  <j>i=2icxJ.  Thus,  by  choosing  <t>2 -<|>1=n/2+rm  or  x2-x1=  (2/1+1)  A/8 ,  the 
in-phase  and  quadrature  components  of  the  signal  can  be  measured.  For  the 
MA86735  device,  the  detector  spacing  is  selected  to  produce  a  nominal  phase 
difference  of  90°  ±15°.  Examination  of  the  output  signals  revealed  that  the  actual 
detector  phasing  was  about  80°  for  the  particular  unit  used  during  the  High  Resolution 
experiment.  A  phase  error  was  also  incurred  due  to  the  fact  that  the  two  detector 
signals  were  not  sampled  exactly  simultaneously.  These  phase  errors  were  corrected 
and  the  amplitudes  of  the  two  detector  signals  were  equalized  during  processing  of  the 
data. 

The  detector  output  signals  were  fed  into  two  parallel  operational  amplifiers  in 
order  to  remove  the  d.c.  bias  and  amplify  the  fluctuating  part  of  the  signal  by  a  factor 
of  100  (20  dB).  The  antenna  used  was  a  standard  gain  horn  antenna  manufactured  by 
Scientific  Atlanta,  which  has  a  beamwidth  of  approximately  12°  and  a  gain  of  22.5  dB 
at  the  operating  frequency  of  10.525  GHz.  The  unit  was  mounted  so  as  to  produce  a 
vertically  polarized  field  with  an  incidence  angle  of  about  45°.  The  data  was  sampled 
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at  a  rate  of  167  sample  pairs  per  second,  using  a  data  acquisition  system  supplied  by 
the  Woods  Hole  Oceanographic  Institution,  and  recorded  on  the  hard  disk  of  a 
personal  computer. 

The  Doppler  radar  was  mounted  on  the  LADAS  catamaran,  along  with  a  video 
camera,  laser  wave  slope  gauge  and  a  mini-SODAR  for  meteorological  observations 
(E.  Bock,  personal  communication).  The  Doppler  radar  was  mounted  next  to  the 
video  camera  on  a  mast  approximately  4  meters  above  the  water  surface,  as  shown  in 
Figure  14.  This  produced  a  footprint  on  the  surface  having  a  diameter  of  about  1 
meter. 

A  total  of  18  hours  of  data  was  collected  during  the  period  13-24  September, 
1991.  Examples  of  the  data  are  shown  in  Figures  15  and  16.  Figure  15  shows  the 
raw  I&Q  signals  recorded  during  a  3.6  second  time  interval  under  low  wind 
conditions  on  September  16.  The  data  was  processed  using  a  short-time  Fourier 
transform  method  to  produce  two-dimensional  (frequency  versus  time)  plots,  as  shown 
in  Figure  16.  These  plots  indicate  that  under  these  low  wind  conditions  the  signal  is 
dominated  by  an  extremely  coherent,  though  time  varying,  component.  Further 
analysis  indicates  that  this  component  is  probably  due  to  a  sidelobe  of  the  antenna 
gain  pattern.  The  third  sidelobe  of  the  horn  antenna  occurs  at  an  angle  of  about  44° 
from  the  boresight  direction,  and  was  therefore  directed  nearly  toward  the  nadir  for 
the  configuration  in  which  the  radar  was  deployed.  The  theoretical  two-way  antenna 
gain  for  this  sidelobe  is  more  than  40  dB  below  the  main  lobe.  Data  from  the  NRL 
4-frequency  radar  (Valenzuela,  1978)  indicates  that  the  radar  cross  section  is  typically 
25  dB  larger  at  nadir  than  at  45°  incidence,  for  vertical  polarization.  However,  under 
low  winds  it  appears  that  this  difference  can  be  much  larger  and  the  return  from  the 
nadir-looking  sidelobe  can  dominate  the  main  lobe  return.  In  future,  it  is 
recommended  that  a  shield  be  installed  below  the  antenna  to  eliminate  this  sidelobe 
return. 
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Hit  C  to  continue  or  S  to  stop 

Figure  15.  Examples  of  Doppler  Radar  In-Phase  and  Qurdrature  Signals 
Recorded  During  the  High  Resolution  Field  Experiment. 
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This  problem  did  not  appear  to  be  present,  at  least  to  the  same  degree,  on 
September  17  during  a  crossing  of  the  ’rip’  feature  which  was  the  focus  of  much 
attention  during  the  first  High  Resolution  experiment.  A  plot  of  the  backscattered 
power  (in  relative  units)  from  the  Doppler  radar  during  one  of  these  crossings  is 
shown  in  Figure  17.  The  average  backscattered  power  during  the  interval  from 
16:17:15  to  16:19:15  EDT  was  five  times  greater  than  the  average  power  during  the 
interval  from  16:16:30  to  16:16:45.  This  increase  is  somewhat  smaller  than  the 
maximum  change  observed  by  the  NRL  airborne  real-aperture  radar  (F.  Askari, 
personal  communication)  for  the  same  or  a  similar  feature.  The  difference  may  be 
due  to  a  difference  in  the  radar  look  direction,  or  may  reflect  some  contamination  of 
the  Doppler  radar  signal  by  the  sidelobe  return  in  the  low-return  area  just  prior  to  the 
rip  crossing. 

Because  of  problems  with  the  WHOI  laser  slope  gauge,  as  well  as  the 
aforementioned  problems  with  the  Doppler  radar,  the  original  intention  of  using  the 
combination  of  these  data  sets  for  model  verification  purposes  has  been  largely 
unfulfilled.  Instead,  a  collaborative  effort  with  Peter  Smith  of  NRL/Stennis  has  been 
initiated,  utilizing  his  measurements  of  the  wave  spectrum  from  a  buoy  as  it  drifted 
across  the  rip  feature,  and  a  comparison  of  model  predictions  with  SAR  observations 
has  also  been  carried  out  as  described  in  the  following  section. 

4.2.  SAR  DATA  ANALYSIS 

Data  was  collected  with  the  P3  SAR  system  on  six  days  during  the  first  High 
Resolution  field  experiment,  on  September  11,  12,  16,  19,  20  and  21.  Of  these  data 
sets,  the  ones  collected  on  the  16®  appear  to  contain  the  largest  number  of  interesting 
features  in  the  vicinity  of  the  research  vessels,  and  these  have  consequently  received 
the  most  attention.  In  fact,  the  only  digitally  processed  data  available  during  the 
period  of  this  project  was  from  data  set  1  (pass  2)  on  the  16®.  Our  analysis  has 
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17  SEPTEMBER  1991 
16:16:25  -  16:20:25  EST 

RECEIVED  POWER 


TIME  (SEC) 

Figure  17.  Backscattered  Power  Received  by  Doppler  Radar  on  September  17 
1991  During  a  Crossing  of  the  ‘Rip’  Feature. 


?>ERIM 


therefore  concentrated  on  this  data  set. 

On  the  morning  of  September  16,  the  USNS  Bartlett  executed  a  box  pattern 
centered  at  approximately  35°  16’N  and  75°  5’W.  The  box  was  oriented  at  about  45° 
from  North  so  that  the  northeast  and  southwest  edges  of  the  box  intersected  the  edge 
of  the  Gulf  Stream  at  approximately  right  angles.  Two  distinct  changes  in  the  water 
temperature  were  observed  aboard  the  Bartlett  at  10:45  and  10:54  EDT  while  heading 
southeast  into  the  Gulf  Stream,  on  the  first  leg  of  the  box  pattern.  The  chief 
scientist’s  log  indicates  the  following  entry  at  10:54:  "Crossing  scum  line. 

Athwartship  speed  jumped  from  below  1  knot  (where  it  has  been  all  this  time)  to  1.8 
then  2.2  knots.  Perturbation  in  gyro  heading  was  about  4  deg,  indicating  a  current 
front  so  sharp  the  auto  heading  device  couldn’t  respond  fast  enough.  STAR  saw 
stratified  water  before  scum  line,  changing  to  isothermal  water  after  it  "  (G. 
Marmorino,  personal  communication). 

SAR  data  set  1  (pass  2)  was  collected  over  this  area  at  approximately  12:08 
EDT.  A  segment  of  this  pass  was  digitally  processed  by  NAWC,  covering  the  area 
shown  in  Figure  18.  In  order  to  compensate  for  the  time  difference  between  the 
overpass  and  the  ship  observations,  the  locations  of  the  water  parcels  sampled  by  the 
ship  at  10:45  and  10:54  EDT  were  projected  using  the  water  velocities  measured  by 
the  Bartlett’s  ADCP,  i.e.  50.8  cm/s  toward  31.4  °T  at  10:45  and  67.6  cm/s  toward 
42.9  °T  at  10:54.  These  projected  locations  are  also  shown  in  Figure  18,  aiong  with 
the  outline  of  a  subset  of  the  digitally  processed  image  which  was  selected  to  include 
these  points.  The  Lw  and  Xw  band  SAR  images  for  this  subset  are  shown  in  Figure 
19.  Based  on  the  position  information,  which  is  of  course  somewhat  uncertain 
because  of  the  relatively  large  time  difference  between  the  shipboard  and  aircraft 
observations,  we  conclude  that  the  two  bright  lines  running  diagonally  through  the 
image  subset  shown  in  Figure  19  correspond  to  the  two  thermal  gradients  observed 
aboard  the  Bartlett,  with  the  rightmost  line  corresponding  to  the  scum  line  and 
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35°6'N 


Figure  18 


Figure  19.  L-Band  and  X-Band  Images  for  Subset  Indicated  in  Figure  18. 

Horizontal  Line  Through  L-Band  Image  is  the  "Double-Nadir"  Return 
Which  Occurs  at  an  Incidence  Angle  of  60°.  Image  Dimensions  are 
170  x  410  Pixels  or  1377  x  3321  Meters,  Each  Pixel  Representing  the 
Sum  of  3  x  5  Original  One-Look  Pixels. 
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velocity  gradient  noted  in  the  chief  scientist’s  log. 

It  is  interesting  to  note  that  the  orientation  of  these  features  (approximately 
70°T)  is  much  different  from  the  measured  current  direction,  which  suggests  the 
possibility  of  a  strong  converging  as  well  as  shearing  current.  The  current  variations 
in  the  vicinity  of  these  features  which  were  noted  in  the  chief  scientist’s  log  were  not 
clearly  observed  by  the  ADCP,  apparently  because  the  horizontal  and/or  vertical 
resolution  of  the  measurements  was  insufficient  to  resolve  the  features.  However,  an 
apparently  similar  feature  (dubbed  the  ’rip’  feature)  was  observed  25  hours  later  by 
both  the  ships  and  the  NRL  airborne  real -aperture  radar  at  around  35°  23 'N  and  75° 
O’W,  i.e.  about  18  km  north  of  the  features  in  the  SAR  image.  These  features  appear 
to  correspond  with  ’shingles’  protruding  from  the  North  edge  of  the  Gulf  Stream, 
which  are  visible  on  AVHRR  thermal  images  and  are  observed  on  these  images  to 
migrate  Northward  about  20  km  per  day.  The  convergence  velocities  in  these 
features  have  been  estimated  from  various  sources  (including  a  sequence  of  NRL 
RAR  images)  to  be  on  the  order  of  25-40  cm/s. 

The  nearest  available  wind  measurement  was  made  aboard  the  R/V  Oceanus  at 
11:30  EDT  on  September  16,  when  the  Oceanus  was  located  at  35°  16’N  and  75° 

O’W.  This  measurement  indicated  a  wind  speed  of  4.2  m/s  from  161T  (J.  Edson, 
personal  communication).  In  attempt  to  match  the  spectral  shapes  measured  by  Jahne 
and  Riemer  (1990)  more  closely  at  low  wind  speeds,  the  net  source  function  discussed 
in  section  2.1  was  modified  by  increasing  the  dissipation  rate  for  wavenumbers  larger 
than  1  rad/cm.  The  angular  distribution  of  the  growth  rate  was  also  modified  so  as  to 
produce  the  ratio  of  upwind  to  downwind  propagating  waves  inferred  by  Plant  and 
Keller  (1990)  from  L-band  Doppler  radar  measurements.  The  current  field  was 
approximated  by  a  linear  variation  in  both  the  normal  and  tangential  components  of 
the  current  from  -20  cm/s  to  +20  cm/s  over  a  10  meter  transition  region. 

Changes  in  the  wave  spectrum  due  to  the  interaction  of  the  ambient  wave  field 
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with  this  current  were  calculated  using  the  model  described  in  Lyzenga  (1991). 
Changes  in  the  radar  backscatter  at  L-band  and  X-band  were  then  calculated  using  the 
two-scale  model  described  by  Plant  (1986).  These  results  are  shown  as  the  solid  lines 
in  Figures  20  and  21.  The  contribution  to  the  backscatter  from  breaking  waves  was 
also  calculated  using  the  model  described  in  sections  2.5  and  3.3  of  this  report.  This 
contribution  is  indicated  by  the  dashed  line  in  Figure  2 1 .  The  slope  variance  used  for 
this  calculation,  as  obtained  from  the  wave-current  interaction  model,  is  shown  in 
Figure  22. 

The  change  in  backscatter  from  the  two-scale  model  (without  wave  breaking 
effects)  is  much  smaller  at  X-band  than  at  L-band,  as  discussed  in  Lyzenga  (1991). 
When  breaking  effects  are  included,  the  changes  in  backscatter  at  X-band  and  L-band 
are  comparable,  for  this  case,  both  being  on  the  order  of  7-8  dB.  This  appears  to 
agree  qualitatively  with  the  SAR  images  shown  in  Figure  19.  A  quantitative 
comparison  is  difficult  because  the  background  signal,  at  least  for  the  X-band  image, 
appears  to  be  below  the  noise  floor.  Plots  of  the  image  intensity  across  the  feature 
(taken  from  the  right-hand  side  of  Figure  19)  are  shown  in  Figures  23  and  24.  From 
these  plots,  the  ratio  of  the  peak  to  background  signal  can  be  inferred  to  be  at  least  a 
factor  of  6  (7.8  dB)  at  X-band  and  a  factor  of  12  (10.8  dB)  at  L-band.  It  may  also 
be  noted  that  some  weaker  features  can  be  seen  in  the  L-band  image  which  do  not 
appear  to  be  present  in  the  X-band  image.  This  is  explainable  in  the  context  of  the 
present  theory  by  reference  to  the  fact  that  the  breaking  wave  return  is  a  strongly 
nonlinear  function  of  the  slope  variance. 
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Figure  20.  L-Band  Radar  Cross  Section  Variations  Across  ‘Rip’  Current  Feature, 
Calculated  From  Wave-Current  Interaction  Model  Combined  With 
Two- Scale  Radar  Backscatter  Model. 
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Wave  Breaking  RCS 

Figure  21.  X-Band  Radar  Cross  Section  Variations  Across  ‘Rip’  Current  Feature, 
Calculated  From  Two-Scale  Model  (Solid  Line)  and  Breaking  Wave 
Model  (Dashed  Line). 
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Figure  22.  Slope  Variances  Calculated  From  Wave-Current  Interaction  Model  for 
‘Rip’  Feature. 
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5.0  SUMMARY  AND  CONCLUSIONS 

During  this  project,  a  new  analytic  solution  of  the  wave  action  equation  was 
derived,  and  various  formulations  were  investigated  for  the  net  source  function  which 
appears  in  this  equation.  The  present  formulation  includes  source  terms  due  to  the 
Phillips  growth  mechanism  and  the  Miles  or  exponential  growth  mechanism,  and 
dissipation  terms  due  to  viscosity,  surfactant  damping,  and  wave  breaking. 

Coefficients  have  been  selected  so  as  to  yield  an  equilibrium  spectrum  which  is  in 
reasonable  agreement  with  observations.  However,  considerable  uncertainty  exists 
with  respect  to  several  of  these  terms  and  the  present  formulation  can  only  be 
considered  as  provisional. 

One  area  in  which  the  present  formulation  of  the  net  source  function  is  quite 
clearly  lacking  is  in  its  neglect  of  nonlinear  interactions,  which  involve  the  transfer  of 
energy  from  one  region  of  the  spectrum  to  another.  As  a  first  step  in  this  direction,  a 
model  for  the  effects  of  wave  breaking  was  developed,  which  describes  the 
contribution  of  breaking  waves  in  terms  of  the  breaking  fraction,  which  is  a  function 
of  the  slope  variance.  This  model  appears  to  predict  the  wind  speed  dependence  of 
the  surface  curvature  spectrum  at  wavenumbers  between  1  and  8  rad/cm.  The 
electromagnetic  scattering  from  breaking  waves  is  calculated  using  the  Kirchhoff 
approximation,  and  the  results  are  found  to  compare  favorably  with  several 
observations  of  the  radar  backscatter  from  breaking  waves. 

The  combined  model  predicts  that  the  radar  backscatter  at  centimeter 
wavelengths  is  strongly  and  nonlinearly  dependent  on  the  slope  variance  associated 
with  wavenumbers  below  about  1  rad/cm.  In  the  open  ocean,  this  variance  is  linearly 
related  to  the  wind  speed  or  friction  velocity.  However,  in  the  presence  of  variable 
surface  currents,  the  slope  variance  may  also  change  dramatically  due  to  the 
interaction  of  surface  waves  with  these  currents.  Comparison  of  the  predictions  of 
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this  model  with  SAR  observations  during  the  first  High  Resolution  field  experiment 
appear  favorable,  but  of  course  many  more  such  comparisons  are  required  in  order  to 
adequately  test  the  model. 
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The  interaction  of  short  surface  waves  with  shearing  and  converging  currents  is  investigated  by 
means  of  an  approximate  analytical  solution  of  the  wave  action  spectral  transport  equation.  Spectral 
perturbations  of  less  than  lOfr  are  predicted  at  centimeter  wavelengths  for  moderately  strong  ocean 
fronts,  although  much  larger  perturbations  are  expected  at  wavelengths  of  the  order  of  1  m.  These 
results,  when  combined  with  a  simple  Bragg  electromagnetic  scattering  model,  do  not  explain  the  large 
backscatter  variations  observed  at  X  band  and  C  band  in  the  vicinity  of  ocean  fronts.  A  two-scale 
electromagnetic  scat'ering  model  provides  an  effective  coupling  between  the  long-wave  spectral 
perturbations  and  the  radar  backscatter  at  small  incidence  angles,  but  this  mechanism  becomes  less 
effective  at  intermediate  incidence  angles.  Observation  of  large  backscatter  variations  at  these  angles 
may  thus  be  an  indication  of  additional  hydrodynamic  and/or  electromagnetic  scattering  mechanisms 
that  are  not  yet  accounted  for. 


1 .  Introduction 

Ocean  fronts  are  characterized  by  large  horizontal  and 
vertical  gradients  in  the  fluid  velocity  as  well  as  density  (i.e., 
temperature  and/or  salinity).  Surface  water  flows  toward  the 
front  and  is  subducted  at  the  frontal  boundary,  causing  a 
converging  current  at  the  surface.  The  component  of  the 
current  parallel  to  the  front  is  also  frequently  observed  to 
have  an  intense  horizontal  shear,  particularly  in  the  case  of 
larger-sole  ocean  fronts  [Garvine  and  Monk.  1974].  These 
two  types  of  surface  current  variations  are  illustrated  in 
Figure  1 .  The  interaction  of  surface  waves  with  these  cur¬ 
rents  causes  a  change  in  the  surface  roughness,  which  also 
implies  a  change  in  the  radar  reflectivity  of  the  surface.  As  a 
result,  linear  features  can  frequently  be  observed  in  radar 
images  of  the  ocean  where  fronts  are  known  to  occur  [e.g., 
Larson  et  al.,  1976;  Mattie  el  al.,  1980;  Hayes ,  1981; 
Veseckv  and  Stewart ,  1982;  Fu  and  Holt.  1983].  Other 
mechanisms  may  also  be  involved  in  some  cases,  including 
variations  in  wind  stress  associated  with  temperature 
changes  across  the  front  and  damping  of  short  surface  waves 
by  surfactant  materials  accumulated  along  the  front. 

In  this  paper  the  wave-current  interaction  mechanism  is 
explored  as  a  possible  explanation  for  the  appearance  of 
fronts  in  radar  images  of  the  ocean  surface.  An  approximate 
solution  of  the  wave  action  spectral  transport  equation  is 
derived  in  the  following  section.  The  properties  of  this 
solution  are  described  in  section  3,  and  the  implications  for 
radar  imaging  of  ocean  fronts  are  discussed  in  section  4, 


2.  Wave-Current  Interactions 

The  interaction  of  surface  waves  with  a  steady,  one¬ 
dimensional  current  field  (i.e.,  a  current  which  is  a  function 
of  only  one  spatial  variable)  can  be  described  in  terms  of  the 
action  spectral  transport  equation 
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where  cgx  is  the  x  component  of  the  wave  group  velocity;  u 
and  v  are  the  x  and  y  components  of  the  surface  current, 
respectively;  kx  and  kY  are  the  x  and  y  components  of  the 
wave  number;  N  is  the  action  spectral  density;  which  is 
defined  as  the  energy  spectral  density  divided  by  the  wave 
frequency;  and  f$(N)  represents  the  net  source  function  for 
wave  action,  which  implicitly  includes  the  effects  of  air-sea 
interactions,  nonlinear  wave-wave  interactions,  viscous  dis¬ 
sipation,  and  other  dissipation  mechanisms  such  as  wave 
breaking  [Phillips,  1984], 

The  form  of  the  net  source  term  on  the  right-hand  side  of 
(1)  is  quite  uncertain.  In  fact,  it  is  not  clear  that  it  can  safely 
be  assumed  to  be  only  a  local  function  cf  N  since  energy 
may  be  exchanged  among  different  regions  of  the  spectrum 
by  nonlinear  interactions.  However,  this  assumption  has 
been  provisionally  adopted  for  the  sake  of  simplicity,  until 
the  importance  of  these  nonlocal  effects  is  demonstrated  and 
a  satisfactory  method  of  accounting  for  them  is  devised.  The 
existence  of  a  stable  equilibrium  condition  in  the  absence  of 
any  currents  implies  that /S(/Vn)  =  0  and  fs(N0)  <  0  where 
Nq  is  the  equilibrium  spectral  density.  Expanding/5(AO  in  a 
Taylor  series  about  this  point,  we  can  then  write 

fs(N)=  -pr(N-  N0)+  y(N-  N0)2  +  ■  ■  ■  (2) 

where  pr  -  -fsiN0)  and  y  =  r/jUVo)-  If  higher-order 
terms  are  neglected,  and  if  the  condition  /s(0)  =  0  is 
imposed,  this  form  is  equivalent  to  that  used  by  Hughes 
[1978],  with  p,  =  -N 0y  =  P  where  p  is  the  initial  growth 
rate  due  to  wind  input.  If  higher-order  terms  are  included, 
the  relaxation  rate  p ,  is  not  necessarily  equivalent  to  the 
growth  rate  p. 

Using  this  expansion  for  the  right-hand  side,  it  is  conve¬ 
nient  to  rewrite  the  spectral  transport  equation  as 

ef 

(*'«  +  u)  —  +  p,f  =  fi(x.  kt.  ky)  +  e,  (3) 
dx 

where  /  =  N/N0  -  1  is  the  fractional  spectral  perturbation 
or  deviation  from  the  equilibrium  state. 


du 

p(x.  k,.  ky) 

dx 


dx  A  1  UN  o 
'  dx  ];V0  t U, 


(4) 
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where  <  -  u.  When,  in  addition,  the  rate  of  shear 

(or  vorticitvt  dx  dx  is  constant,  this  reduces  further  to 


M  I 

/‘-Mi '  -  —if 


Fig.  la 


dtt  d\  ~  U.  dx  dx  -  constant 

For  the  case  in  which  the  strain  rate  du  dx  is  constant  but 
not  zero.  (3)  no  longer  has  constant  coefficients,  but  a 
solution  can  still  be  found  by  using  the  integrating  factor 

AUi  =  exp  {0‘  In  (c,,  -  i/lf  =  (c9,  -  // 1 e  (8) 

where  0’  =  fSrldnidx)  ~ 1 .  Assuming  that  the  shear  rate 
dvidx  is  also  constant,  the  solution  can  be  written  tn  the 
form 


M  A(x0) 


0rj  A  t  jr ) 


duldx  =  constant,  dvidx  =  constant. 


ft  is  also  instructive  to  consider  two  limiting  cases  relative 
to  the  relaxation  rate  0,.  In  the  "lossless"  limit,  i.e..  for 
0,  =  0.  the  above  solution  for  du'dx  *  0  approaches 


Fig.  16 

Fig.  I .  Illustration  of  surface  current  variations  in  the  vicinity  of 
an  ocean  front:  (a)  converging  current  [duldx  <  0.  dvidx  =  0)  and 
lb)  shearing  current  (duldx  =  0.  dvidx  <  0). 


fix)  =/U0)  +  — 


which  reduces,  for  the  case  duldx 


0,  =  0.  duidx*  0 


fix)  =  /(xn)  +  jut 


0,  =  0.  du  dx  =  0. 


M=  U , 


dx  )  '  A'n  bkx 


yKnf~ 


higher-ordet  terms.  (5) 


fix  1  = 


H 


fXt  X.  k  ,  ,  f,|(’ 


In  the  opposite  limit,  as  0r  ■ 
approach 


both  of  the  above  solutions 


If  the  current  gradients  duldx  and  dvidx  are  small,  the 
spectral  perturbation / is  expected  to  be  small  as  well.  Thus 
tr  contains  only  second-order  and  higher-order  terms.  We 
can  therefore  hope  to  get  some  insight  into  the  problem  by 
considering  the  solution  of  the  equation  in  which  r,  is 
neglected.  The  validity  of  this  approximation  in  some  cases 
of  interest  when  the  spectral  perturbation  is  not  small  is 
discussed  in  the  following  section. 

Even  when  r,  is  neglected,  (3)  is  difficult  to  solve  in  the 
general  case  because  of  the  nonconstant  coefficient  in  the 
first  term.  However,  there  are  several  special  cases  for 
which  exact  solutions  of  the  linearized  equation  exist.  One  of 
these  ts  the  ca-e  when  u  is  constant,  i.e..  the  strain  rate  duldx 
is  zero,  so  the  only  current  gradients  are  in  the  flow  pantile!  to 
the  front,  as  illustrated  in  Figure  16.  In  this  case.  (3)  docs  have 
constant  coefficients  and  can  be  readily  solved  by  using  the 
integrating  factor  exp  {-f},xHcgr  +  «)}.  to  yield 


fix)  =  *z'0r  0, » (c9,  +  uML  <l2t 

where  L  is  a  characteristic  length  scale  for  the  current 
pattern.  This  solution  is  obvious  by  inspection  of  (3).  and 
furthermore,  it  can  be  seen  to  hold  also  at  the  locations 
where  rflI  *  u  -  0  for  any  nonzero  0r.  This  is  the 
well-known  "blocking”  or  “resonance"  condition. 

Although  derived  for  a  constant  or  uniform  current  gradi¬ 
ent,  the  above  solutions  can  be  applied  to  an  arbitrary 
current  field  by  breaking  it  up  into  a  set  of  piecewise  linear 
segments.  For  the  purpose  of  exploring  the  general  behavior 
of  the  wave-current  interactions  near  an  ocean  front,  how¬ 
ever.  the  current  field  can  be  approximated  by  the  simple 
ramp  functions 


II  i  it  -  -u,.\i  m  -n  <  a  <  >v 


u i  a  i  =  -II, 


i  t  X  I  =  -  ;  ,V  U  -  IV  <  A  <  H 


du’dx  =  0  (6l 


LVZi  NGA:  INTERACTION  Ot  WAVES  WlTH  CX'tAN  h RUNTS 


10. Th’ 


which  imply  a  strain  rate  dmdx  =  -u,  <»  and  a  shear  rate 
dvidx  ~  -vjw  in  the  central  region. 

The  solution  is  evaluated  at  any  given  location  by  first 
determining  the  appropriate  value  of  a,,  on  the  basis  of  the 
sign  of  c9<  -r  it.  For  cgx  -  u  >  0  the  boundary  condition  is 
applied  at  the  left  boundary  of  the  region  under  consider¬ 
ation,  i.e..  at  a o  =  w  for  points  inside  the  region 
-w  <  x  <  k  and  at  a0  =  h  for  points  in  the  region  a  >  «  . 
The  appropriate  boundary  value  for  this  case  is  fix )  =  0  for 
as  -  u\  For  cgx  +  u  <  0  the  boundary  condition  is  applied 
at  the  right  boundary,  and  the  boundary  value  is  fix )  =  0  for 
x  a  tv. 


Suik.  <b  I  —  f  ,|td>iA  (20l 

where  p  =  4  in  the  gravity  region  and  p  ~  8  in  the  capillary 
region.  The  angular  dependence  of  the  spectrum  may  be 
represented  by  the  function 

,  /  <6  -  <6  „  \ 

F„ut>)  =  <„  cos-”  I — - - j  121) 

where  n  ~  2-5  in  the  intermediate  region  of  the  spectrum, 
and  n  =  10  near  the  peak  of  the  spectrum  [Pierson  and 
Stacy.  1973].  Using  these  forms,  the  factor  p  in  (4)  can  be 
written  as 


3.  Spectral  Perturbations  Induced 
by  Ocean  Fronts 

The  spectral  perturbations  caused  by  the  shearing  and 
converging  currents  in  an  ocean  front  depend  on  the  wind 
speed  and  direction  (through  the  relaxation  rate  /3r  and  the 
equilibrium  spectrum  ;V0)  as  well  as  the  current  field.  Using 
the  Hughes  [1978]  formulation  for  the  net  source  function, 
the  relaxation  rate  can  be  considered  to  be  approximately 
equivalent  to  the  initial  growth  rate  0.  On  the  basis  of  wave 
tank  measurements  as  well  as  theoretical  considerations. 
Plant  and  Wright  [1977]  proposed  a  growth  rate  of  the  form 

0  =  0.0Aiu,/c)2w  (15) 

where  w*  is  the  friction  velocity  (which  is  roughly  1/30  of  the 
wind  speed  U  measured  at  a  standard  height  of  19.5  m  above 
the  surface),  c  is  the  phase  velocity,  and  at  is  the  frequency 
of  the  waves.  The  dependence  of  /3on  the  wave  propagation 
direction  relative  to  the  wind  direction  could  not  be  observed 
in  the  linear  wave  tank  used  for  these  measurements.  How¬ 
ever,  Plant  [1986]  has  suggested  an  angular  dependence  of 
the  form 


Fi<t>)  =  cos  (<6  -  d>H.)  | <t>  -  <6»j  <  w/2 

F(<6)  =  0  | <t>  -  <6* I  >  ir/2 


(16) 


where  <6„  is  the  wind  direction 
The  equilibrium  spectrum  enters  this  solution  through  its 
logarithmic  gradient 


k  6Nq  k  1  dA/Q 

A'o  dkx  ~  N0  dk  C°S  *  N0  dd> 


(17) 


where  /V0  is  the  action  spectral  density,  which  is  related  to 
the  energy  spectrum  E0( A ,  <t>)  and  the  wave  height  spectrum 
Soik.  <t>)  through  the  equation  [ Phillips ,  1980} 


1  Cl> 

iV„( k.  <6>  =  -  E0(k.  <61  =  p-S0(k,  <6)  (18) 

at  k 


where  k  is  the  magnitude  of  the  wave  number  and  at  is  the 
intrinsic  wave  frequency,  which  is  given  by 

at2  =  gk  +  iTlp)kX  (19) 

where  g  is  the  gravitational  acceleration,  T  is  the  surface 
tension,  and  pis  the  density  of  water  (T/p  =  74  cm’/s~  for 
a  clean  water  surface). 

For  wave  numbers  much  larger  than  the  peak  wave 
number  k,t  -■  g KJ2 .  the  equilibrium  wave  height  spectrum 
can  be  modeled  as  a  power  law,  i.e.. 


/  du  dv 

p  =  —  cos  <6  -  —  sin  <6 
\dx  dx 


-( p  +  1  -  Cglc)  cos  <6 


+  n  tan 


sin  <6J 


(22) 


The  factor  p.  is  therefore  of  the  same  order  of  magnitude  as 
the  current  gradients,  which  are  in  the  range  of  0.01-0.001 
s-1  for  moderately  strong  ocean  fronts.  The  relaxation  rate 
is  typically  of  the  order  of  0. 1-1.0  s'1  for  wavelengths  from 
1  to  10  cm  and  moderate  wind  speeds.  Thus  the  fractional 
spectral  perturbations  predicted  by  this  model  are  typically 
only  a  few  percent  in  this  wavelength  range.  On  the  other 
hand,  much  larger  fractional  perturbations  are  predicted  for 
longer  wavelengths,  where  the  relaxation  rate  is  much 
smaller. 

At  this  point,  before  presenting  some  examples  to  illus¬ 
trate  the  behavior  described  in  the  previous  paragraph,  it  is 
worth  mentioning  two  problems  that  occur  when  the  relax¬ 
ation  rate  (3,  is  equated  with  the  growth  rate  /3  as  discussed 
at  the  beginning  of  this  section.  The  first  problem  is  associ¬ 
ated  with  the  angular  dependence  of  /3.  Using  the  angular 
dependence  shown  in  (16),  the  growth  rate  falls  to  zero  for 
waves  traveling  at  angles  larger  than  90°  with  respect  to  the 
wind  direction.  If  the  same  angular  dependence  is  assumed 
for  the  relaxation  rate  /3r,  unrealistically  large  spectral 
perturbations  are  predicted  at  these  angles,  as  shown  by 
( I0H 12).  The  other  problem  occurs  when  the  growth  rate  is 
corrected  for  viscous  dissipation  effects.  The  proper  correc¬ 
tion  to  the  growth  rate  is  to  subtract  4vk‘  from  (3.  where  vis 
the  kinematic  viscosity  [ Plant  and  Wright.  1977],  However, 
if  this  is  done,  the  net  growth  rate  becomes  negative  at  high 
wave  numbers  or  at  large  angles  to  the  wind.  Obviously,  the 
relaxation  rate  cannot  be  simply  equated  to  the  net  growth 
rate  in  these  spectral  regions,  since  the  relaxation  rate  must 
be  positive. 

The  solution  to  both  of  these  problems  is  to  use  a  more 
consistent  model  for  the  net  source  function  and  to  calculate 
the  equilibrium  spectrum  by  solving  the  equation /5(,Vfl)  = 
0.  Steps  in  this  direction  have  been  taken,  for  example,  by 
Donelan  and  Pierson  [1987]  and  Plant  [1986],  The  simplest 
such  model  would  consist  of  a  wind-forcing  term  propor¬ 
tional  to  \  and  a  dissipation  term  proportional  to  Nr  where 
p  >  1 .  If  viscous  dissipation  is  included  in  such  a  model,  the 
relaxation  rate  (ir  -  f si  A/0)  still  vanishes  as  the  net  growth 
rate  ifi  -  4i ik2)  approaches  zero,  leading  to  large  fractional 
perturbations.  However,  the  equilibrium  spectrum  also  ap¬ 
proaches  zero  at  these  points,  so  the  perturbed  spectrum 
i V  =  /V0<  1  +  /)  remains  finite,  at  least  if  p  <  2. 
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A  somewhat  more  satisfactory  result  is  cbtained  if  another 
source  term  which  is  independent  of  S'  is  added  to  f\ ( ,V ) .  In 
that  case  the  equilibrium  spectrum  and  the  relaxation  rate 
are  both  positive  for  all  wave  numbers.  The  equilibrium 
spectrum  is  given  by  the  solution  of  the  equation 

/s< AV  =  U3  -  -  yAtJ  +  II  =  0  (23) 

where  II  is  the  constant  term  (which  physically  corresponds 
to  the  Phillips  resonant  growth  mechanism),  and  the  relax¬ 
ation  rate  is  then  given  by 

,  n 

0,  =  ~fS(N o)  =  (P  “  1K0  ~  4vA‘)  +  p—.  (24) 

tVo 

For  moderate  wind  speeds,  wavelengths  larger  than  a  few 
centimeters,  and  wave  propagation  directions  within  90°  of 
the  wind  direction,  the  constant  term  in  (23)  is  negligible  in 
comparison  with  the  first  two  terms,  and  it  follows  that 

P,~{p-  l)(/3  -4vk2)~(p-  l)/3  (25) 

in  this  region.  On  the  other  hand,  when  4 uk2  >  /3,  the 
equilibrium  state  is  determined  by  the  balance  between  the 
first  and  last  terms  in  (23),  so  that  N0  =  n/(4j-A2  -  )3)  and 

pr  =  4vk2-  p  ~4vk2.  (26) 

The  relaxation  rate  thus  appears  to  be  the  larger  of  ( p  -  1  )0 
and  4  vk  - ,  at  least  in  the  extreme  cases  (what  happens  in  the 
intermediate  region,  where  (p  -  l)/3  *  4vk2,  is  not  yet 
clear).  For  order-of-magnitude  estimates  we  have  therefore 
approximated  the  relaxation  rate  as 

/3,  =  (p-  1)0  +  4vA2  (27) 

with  p  =  2.  To  avoid  the  discontinuity  in  the  derivative  of 
)3(d>)  at  <b  -  =  rr/2  implied  by  (16)  and  to  allow  for  some 

variability  in  the  wind  direction,  we  have  also  chosen  the 
angular  dependence 

,  U  ~  <t>„\ 

Fld>)  =  cos'  I - - - I  (28) 

for  the  growth  rate.  Finally,  in  lieu  of  a  complete,  self- 
consistent  formulation  of  the  net  source  function  as  dis¬ 
cussed  above,  we  have  used  the  equilibrium  spectrum  pro¬ 
posed  by  Bjerkaas  and  Riedel  [  1979]  with  a  cos4(d>/2)  angular 
dependence. 

To  illustrate  the  general  behavior  of  the  solution,  some 
results  were  generated  using  the  current  pattern  specified  by 
(13)  and  (14)  with  uc  =  r,  =  5  cm/s  and  *»•  =  10  m.  for  a 
wind  speed  of  5  m/s  along  the  y  axis,  parallel  to  the  front.  A 
contour  plot  showing  the  wave  number  dependence  of  the 
spectral  perturbation  at  the  center  of  the  current  pattern 
(x  =  0)  is  presented  in  Figure  2.  The  radial  distance  on  this 
plot  represents  the  logarithm  of  the  wave  number,  and  the 
angle  corresponds  to  the  W'ave  propagation  direction.  Thus 
points  falling  along  a  vertical  line  through  the  center  of  the 
plot  represent  waves  traveling  parallel  to  the  front.  The 
dashed  curves  represent  the  locus  of  points  with  a  given 
wavelength  as  indicated  on  the  plot.  The  solid  curves  repre¬ 
sent  linearly  spaced  contours  of  the  dimensionless  quantity 

A4f.S(A.  <*)  -  Snik.  <M]  =  k*S0(k.  4>)/U .  (29) 


Fig.  2.  Contour  plot  of  the  spectra)  perturbation  function  for 
duldx  =  dvldx  =  -0.005  s " 1  with  a  wind  speed  of  5  m/s  parallel 
to  the  front. 


where  50  is  the  equilibrium  height  spectrum  and  /  is  the 
fractional  spectral  perturbation.  The  peak  value  of  this 
function  occurs  at  a  wavelength  of  approximately  1  m  and  a 
wave  propagation  direction  of  about  60°  relative  to  the  x 
axis.  A  minimum  occurs  at  nearly  the  same  wavelength  for  a 
wave  propagation  direction  of  about  108°,  and  another 
maximum  occurs  for  slightly  shorter  waves  propagating 
almost  directly  into  the  front,  from  right  to  left.  Plots  of  the 
fractional  spectral  perturbations  at  these  wavelengths  and 
directions  versus  position  are  shown  in  Figure  3. 

The  wave  number  dependence  of  the  spectral  perturbation 
function  shown  in  Figure  2  is  related  to  the  wave  number 
dependence  of  N0  and  fir  and  is  also  a  result  of  the  combined 
effects  of  the  shearing  and  converging  currents  in  the  front. 
If  there  were  only  converging  currents,  as  shown  in  Figure 
la.  the  spectral  perturbation  would  be  an  even  function  of 
A,  for  a  wind  direction  parallel  to  the  front.  On  the  other 
hand,  if  there  were  only  shearing  currents,  as  shown  in 
Figure  16,  the  spectral  perturbation  would  be  an  odd  func¬ 
tion  of  kx,  with  positive  values  for  k,  >  0  and  negative 
values  for  kx  <  0  for  the  geometry  assumed  in  this  example. 
Thus  the  large  peak  at  <t>  =  60°  in  Figure  1  is  due  to  a 
reinforcement  of  the  effects  of  shear  and  convergence,  while 
the  smaller  peak  at  <S>  -  180°  and  the  minimum  at  <t>  =  108° 
are  due  to  a  partial  cancellation  of  these  effects. 

Figure  3  shows  that  the  fractional  spectral  perturbation  is 
larger  than  I  at  some  locations,  which  raises  a  question  as  to 
the  validity  of  the  linearized  solution  used  for  these  calcula¬ 
tions.  In  order  to  address  this  question,  the  “source"  term 
fiix.  k,,  ky)  and  the  “error"  term  Er(x,  k,.  ky )  in  (3)  were 
calculated  and  are  plotted  in  Figures  4-6  for  the  same 
wavelengths  and  directions  as  used  in  Figure  3.  Figure  4 
shows  that  er  becomes  an  appreciable  fraction  of  n  as  the 
spectral  perturbation  approaches  I.  However,  since  the 
spectral  perturbation  depends  on  the  integral  of  these  terms, 
the  accumulated  effect  of  the  error  term  is  only  about  2 59c  at 
the  nght-hand  boundary  of  the  current  pattern.  Note  that  f, 
and  (i  have  the  same  sign  within  the  current  pattern,  so  the 
linearized  equation  underestimates  the  spectral  perturbation 
in  this  region.  Outside  the  region  where  the  current  is 
varying  (i.e.,  for  x  >  10  m).  er  is  negative,  which  implies 
that  the  full  solution  decays  somewhat  more  rapidly  than  the 
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linearized  solution  in  this  region.  These  conclusions  have 
been  confirmed  by  spot  comparisons  with  "exact"  calcula¬ 
tions  using  the  numerical  model  described  by  Lyzenga  and 
Bennett  [1988]. 

A  similar  analysis  of  Figures  5  and  6  shows  that  the 
magnitude  of  the  negative  perturbations  in  the  vicinity  of 
<t>  =  108°  are  slightly  overestimated,  and  the  positive  pertur¬ 
bations  around  <t>  =  180°  are  slightly  underestimated  by  the 
linearized  model.  For  wave  numbers  well  away  from  these 
peaks,  e,  is  typically  much  smaller  than  p,  indicating  that  the 
linearized  model  performs  quite  well,  as  expected.  The 
spectral  perturbations  predicted  by  this  model  at  a  wave¬ 
length  of  5  cm  are  well  under  10%  for  all  ang(es  of  propaga¬ 
tion. 


4.  Radar  Imaging  of  Ocean  Fronts 


The  scattering  of  microwave  radiation  from  the  ocean 
surface  is  commonly  assumed  to  be  described  by  the  Bragg 
scattering  model,  according  to  which  the  backscattering 
radar  cross  section  per  unit  area  is  given  by 

<7O(0,  4>)  =  8ir*$G(0)r(0)  (30) 

where  8  is  the  incidence  angle  and  4>  is  the  azimuthal  look 
angle,  k0  is  the  electromagnetic  wave  number,  G(0)  is  a 
polarization-dependent  geometric  factor  given  by 


G„(0)  = 


cos'*  6 


Fig.  4.  Magnitude  of  the  source  term  p(  x.  k,.  ky)  and  error  term 
e,(jr,  k„.  kv)  for  A  =  0.99  m  and  <t>  =  60°. 


[1  +  (1/ye)  COS  0]4 
for  horizontal  polarization  and 


SPECTRAL.  PERTURBATION 


Gv(0>  = 


cos4  0(1  +  sin’  9 )' 
[cos  8  +  (I/\c)]4 


for  vertical  polarization,  where  e  »  \  is  the  relative  dielec¬ 
tric  constant  of  seawater,  and 

f(0)  =  S(kB,  <t>)  +  S(kB.  <t>  +  rr)  (33) 


Fig  3  Fractional  spectral  perturbation  versus  position,  at  the 
three  wavelengths  and  directions  indicated  in  legend. 


where  S(k ,  <f>)  is  the  surface  elevation  spectrum  and  kB  = 
2ko  sin  9  is  the  Bragg  resonant  surface  wave  number. 

Combining  this  model  directly  with  the  results  described 
in  the  previous  section,  we  would  expect  ocean  fronts  to  be 
very  rarely  detectable  by  radars  operating  at  wavelengths  in 
the  I-  to  10-cm  range,  since  the  predicted  changes  in  the 
surface  spectrum  at  the  Bragg  wavelength  are  typically  only 
a  few  percent.  However,  several  instances  of  such  detection 
have  been  reported  [ Johannessen  el  al.,  1991;  F.  Askari  et 
al..  An  estuarine  front  viewed  by  an  imaging  radar,  submit¬ 
ted  to  Journal  of  Geophysical  Research.  1991]  for  imaging 
radars  operating  at  C  band  (5  cm)  and  X  band  (3  cm).  The 
fractional  changes  in  backscatter  for  these  cases  are  of  the 
order  of  1  or  larger.  If  the  wave-current  interaction  model 
described  above  is  assumed  to  yield  an  adequate  represen¬ 
tation  of  the  changes  in  surface  roughness  associated  with 
these  fronts,  then  the  dominant  effects  on  the  radar  back¬ 
scatter  must  be  due  to  surface  waves  longer  than  the  Bragg 
wavelength. 

The  effects  of  longer  surface  waves  can  be  included  in  the 
Bragg  model  by  means  of  a  procedure  described  heuristically 
by  Wright  [1968]  and  rederived  by  Brown  (1978).  Valenzuela 
[1978],  Thompson  [1988],  and  others,  using  various  methods. 
Intuitively,  these  effects  can  be  thought  of  as  being  due  to 
changes  in  the  local  angle  of  incidence  caused  by  the  tilting 
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i  / <i<70\  c’( Hi  r'( 0 ) 

m  ( 0 )  =  —  -  = - -  - 

o0\  He  I  Gun  Ft  6*  > 

as 


(37) 


1  # 

Ml  = - [mar0]  =  m'(8 )  ~  m-(0  ).  (38) 

O' o  00 

Using  the  above  expressions  for  Gw<0)  and  G\(8),  the 
polarization-dependent  part  of  m(8)  can  be  written  as 

G'n(6i  4  sin  8 

mH(  8 )  =  — - =  -= - 4  tan  8  (39) 

Gf/iO)  \e  +  cos  8 

for  horizontal  polarization  and 

Gv(0)  4  sin  6  8  sin2  8  tan  8 

mv(0)  =  — — —  - = - - -  (40) 

Gv{0)  (l/Nf)  +  cos  9  1  +  sin*  6 

for  vertical  polarization.  Finally,  if  the  symmetrized  spec¬ 
trum  can  be  expressed  as  a  power  law  in  k ,  i.e.„ 


<ro{0.  4>)  =  vo(0,  <t>) 
where 

I'd -Vo 


I  + 


(vf) 


Pi 


(v:) 


1 


Ml  = 


Col  <30- 


G”(6) 

G(8) 


G'(8)  r<«)  nei 
+  2  - - - + 


G{e)r(6)  ne> 


and 


M; 


1  fd-<70\  2(2\/?-l) 


(34) 


(35) 


(36) 


<r0\  0a  "  /  sin*  0 

where  a  =  tan'1  tj2  is  the  tilt  angle  perpendicular  to  the 
plane  of  incidence,  R  =  GV!GH  for  horizontal  polarization, 
and  R  =  GH!GV  for  vertical  polarization  [Plant,  1988). 

The  sensitivity  factor  mi  can  be  expressed  in  terms  of  the 
tilt  modulation  transfer  function 


S(k,  <{>)  +  S(/t,  <t>  +  it)  =  M4>)k'p 


(41) 


i_Ev^cND 
soiree 
error  te^r 

Fig.  5.  Magnitude  of  the  source  term  x,  kx ,  ky)  and  error  term 
Br{x t  kx.  ky)  for  A  =  0.99  m  and  4>  =  108°. 


of  the  surface  by  the  longer  waves.  Since  the  radar  cross 
section  is  a  nonlinear  function  of  the  incidence  angle,  aver¬ 
aging  over  all  the  slopes  occurring  on  the  surface  results  in  a 
net  change  in  the  mean  backscattered  power,  even  though 
the  mean  slope  is  zero. 

Using  this  approach,  an  expression  for  the  radar  cross 
section  can  be  derived  in  terms  of  the  long-wave  slope 
variances  (p2 )  and  <tj22)  in  the  plane  of  incidence  and  in  the 
perpendicular  direction,  respectively  [Plant,  1986].  A 
change  in  r)(  is  equivalent  to  a  change  in  the  local  incidence 
angle  6,  while  a  change  in  rj:  causes  a  change  in  the 
polarization  of  the  incident  radiation  relative  to  the  local 
normal.  Assuming  the  mean  slope  is  zero,  the  mean  radar 
cross  section  including  these  effects  to  second  order  in  the 
surface  slope  is 


then  the  polarization-independent  part  of  m{0)  can  be  writ¬ 
ten  as 


ms{0)  =  r'(0)/r(<?)  =  -p  cot  0. 


(42) 


Perturbations  in  the  wave  spectrum  can  thus  be  seen  to 
affect  the  radar  backscatter  in  three  ways,  according  to  this 
model:  (1)  through  changes  in  the  spectral  density  at  the 
Bragg  wave  number.  (2)  through  changes  in  the  large-scale 
slope  variances  <tj2>  and  (rj|>.  and  (3)  through  changes  in  the 


Fig.  6. 


Magnitude  of  the  source  term  m<t.  k, .  k, )  and  error  term 
kg,  k ,)  for  A  =  0.55  m  and  <b  =  180° 
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gradient  of  the  spectrum  near  the  Bragg  wave  number  li  e., 
the  value  of  p). 

On  the  basis  of  the  predictions  of  the  wave-current  interac¬ 
tion  model  discussed  in  section  3.  the  second  of  these  would 
appear  to  be  the  dominant  mechanism  for  imaging  ocean  fronts 
at  the  higher  microwave  frequencies  (C  and  X  band),  since  only 
small  changes  in  the  wave  spectrum  are  predicted  at  the  Bragg 
wave  numbers  for  these  radar  frequencies.  The  values  of  fj.,12 
and  are  typically  of  the  order  of  10-100,  depending  on  the 
incidence  angle  and  polarization  [Plant,  1986,  1988],  Thus 
changes  in  the  slope  variance  of  the  order  of  0.01  can  cause 
substantial  changes  in  the  radar  cross  section  and  could  explain 
at  least  qualitatively  the  appearance  of  frontal-related  features 
in  such  images. 

For  the  example  case  discussed  in  section  3.  the  maximum 
change  in  the  X  band  radar  cross  section  at  a  20°  incidence 
angle  is  about  1  dB,  as  shown  in  Figure  7.  Such  a  feature 
would  be  visible  within  the  speckle  or  Rayleigh  noise  back¬ 
ground  of  a  radar  image.  However,  the  magnitude  of  this 
change  decreases  rapidly  with  increasing  incidence  angle, 
particularly  for  vertical  polarization,  as  shown  in  Figures  8 
and  9.  The  backscatter  variations  shown  in  these  plots  would 
be  barely  visible  or  not  detectable,  depending  on  the  radar 
resolution.  On  the  other  hand,  the  backscatter  modulations 
observed  in  the  existing  data  sets  appear  to  be  roughly 
independent  of  incidence  angle.  Consequently,  it  is  not  clear 
that  the  type  of  model  described  in  this  paper  is  capable  of 
adequately  explaining  the  observed  backscatter  modulations 
at  intermediate  incidence  angles,  even  though  it  appears  to 
do  so  at  small  incidence  angles. 

5.  Conclusions 

Changes  in  surface  roughness  can  be  caused  by  the 
interaction  of  surface  waves  with  both  the  converging  and 

'  crp*  '■y,- 
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Fig  7  Backscatter  variation  predicted  by  two-scale  model  at  X 
band.  20°  incidence  for  test  case  discussed  in  text. 
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Fig.  8.  Backscatter  variation  predicted  by  two-scale  model  at  X 
band,  30°  incidence  for  test  case  discussed  in  text. 


the  shearing  currents  associated  with  ocean  fronts.  The 
largest  spectral  perturbations  typically  occur  at  wavelengths 
of  the  order  of  1  m.  The  perturbations  decrease  at  longer 
wavelengths  when  the  group  velocity  becomes  much  larger 


Fig  9.  Backscatter  variation  predicted  by  two-scale  model  at  X 
band.  40°  incidence  for  test  case  discussed  in  text. 
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than  the  current  speed  and  at  shorter  wavelengths  where 
relaxation  effects  play  a  predominant  role.  The  shape  of  the 
spectral  perturbation  function  depends  on  the  wind  direction 
and  is  also  influenced  by  the  nature  of  the  currents,  with  the 
effects  of  shearing  and  converging  currents  tending  to  rein¬ 
force  each  other  in  certain  spectral  regions  and  to  cancel 
each  other  in  other  regions. 

Because  of  the  small  spectral  perturbations  predicted  at 
centimeter  wavelengths,  simple  Bragg  scattering  does  not 
appear  to  be  sufficient  to  explain  the  appearance  of  fronts  in 
radar  images  collected  at  C  band  and  higher  frequencies.  A 
qualitative  explanation  is  provided  by  the  two-scale  scatter¬ 
ing  model,  which  introduces  a  coupling  between  the  long¬ 
wave  spectral  perturbations  and  the  radar  backscatter.  How¬ 
ever.  this  coupling  is  relatively  weak  at  intermediate 
incidence  angles,  and  it  is  not  yet  clear  that  the  observed 
backscatter  variations  at  these  angles  are  adequately  ex¬ 
plained  by  this  model.  Consequently,  additional  hydrody¬ 
namic  and/or  electromagnetic  effects  which  are  not  presently 
accounted  for  may  be  important  in  some  cases.  For  example, 
nonlinear  wave-wave  interactions  or  wave  breaking  may 
result  in  the  cascade  of  energy  from  the  longer  waves  down 
to  centimeter  wavelengths,  necessitating  additional  nonlocal 
source  terms  in  the  wave  action  equation.  Furthermore,  the 
surface  statistics  associated  with  highly  nonlinear  or  break¬ 
ing  waves  may  violate  the  assumptions  of  the  two-scale 
model  and  require  a  reexamination  of  the  scattering  prob¬ 
lem. 
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